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FOREWORD 


This  Project,  (6906)  "Nuclear  Weapon  Effects  on  Space  Vehicles"  and  Task,  (690601)  "Deter¬ 
mination  of  High  Altitude  Nuclear  Weapon  Effects  on  Space  Vehicles,  "  are  a  part  of  the  Air  Force 
Systems  Command  applied  research  program  710A,"  NUCLEAR  WEAPONS  EFFECTS.  " 

The  study  was  initiated  by  the  Flight  Dynamics  Laboratory,  WWRMD,  under  Project  No.  6906, 
Task  No.  690601.  The  project  officer  was  Mr.  L.  E.  Gilbert.  The  investigation  was  conducted  by 
the  Southwest  Research  Institute  during  the  period  from  1  March  1961  through  31  March  1962. 
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ABSTRACT 


Procedures  for  analytically  predicting  the  response  of  missile  bodies  to  blast  loadings  are 
presented.  The  investigation  involves  the  behavior  of  cylindrical  shells  (with  various  end- closures) 
and  circular,  flat  plates.  The  numerical  results  obtained  from  the  analytical  methods  compare 
favorably  with  the  experimental  data  acquired  during  the  study. 
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1.  INTRODUCTION 


In  evaluating  the  vulnerability  of  missile  systems  to  nuclear  antimissile  weapons,  one 
aspect  of  considerable  importance  is  the  response  of  the  vehicle's  structure.  The  analytical  (and 
to  a  lesser  extent,  the  experimental)  information  presented  in  this  report  is  directed  towards 
developing  procedures  for  predicting  the  response  of  missile  bodies  as  represented  by  cylindrical 
shells  with  various  types  of  end- closures . 

Within  the  concept  of  structural  response,  the  most  significant  parameters  which  con¬ 
stitute  a  measurement  of  the  damage  a  missile  will  sustain  as  the  result  of  blast  loadings  are  the 
deformations  and  accelerations.  Consequently,  each  of  the  analytical  procedures  set  out  to  ini¬ 
tially  define  the  time  dependent  variation  of  the  displacement  components.  Once  this  is  accom¬ 
plished,  it  becomes  a  simple  matter  to  obtain  accelerations  as  functions  of  time.  For  the  particular 
shock  sensitive  components  associated  with  each  missile  system,  this  would  provide  an  indication  of 
the  critical  acceleration  forces  to  be  expected. 

It  is  analytically  expeditious  to  subdivide  the  overall  missile  structure  into  two  parts,  each 
representing  different  vulnerability  aspects.  The  first  is  to  consider  the  cylinders'  end-closures 
only  insofar  as  their  restraint  contributes  to  deformational  response  of  the  shell.  Here,  it  is 
implied  that  the  missile  shell  itself  is  most  vulnerable  to  nuclear  weapon  effects.  The  other 
approach  is  to  consider  the  cylinder  end-closures  themselves.  For  example,  a  flat,  circular  plate 
would  require  careful  analysis  if  it  can  be  shown  its  response  is  considerably  more  severe  than  that 
experienced  by  the  shell. 

Accordingly,  in  Section  2,  comparatively  straightforward  analytical  procedures  for 
acquiring  displacements  that  are  reasonably  accurate  and  not  unnecessarily  laborious  are  pre¬ 
sented.  For  the  cylindrical  shell,  both  the  small  and  large  deflection,  linear-elastic  theories  are 
used.  In  the  small  deflection  approach,  the  effect  of  end-closures  {simple,  fixed  and  elastic  sup¬ 
port)  are  considered.  Two  treatments  of  the  circular,  flat  plate  are  presented:  one  entails  the 
dynamic  response  of  circular  plates  at  large  amplitudes  and  the  other  the  plastic  collapse  of  cir¬ 
cular  plates  under  blast  loadings.  Each  analytical  procedure  is  carried  out  in  sufficient  detail  so 
as  to  present  results  which  are  immediately  applicable  to  specific  problems. 

In  Section  3,  extensive  numerical  computations  are  presented  for  the  analytical  procedure 
based  on  the  linear- elastic ,  small  displacement  shell  theory.  In  evaluating  the  analytical  pro¬ 
cedures  for  predicting  the  response  of  cylindrical  shell  structures  to  blast  loadings,  there  are 
several  logical,  as  well  as  practical,  reasons  for  initially  investigating  the  adequacy  of  the  simplest 
method.  Foremost  among  these  reasons  is  the  fact  that  these  methods  provide  numerical  solutions 
for  a  wide  variety  of  variables  with  a  minimum  expenditure  in  computational  effort.  The  extent  to 
which  various  load  and  time  parameters  and  certain  approximations  influence  the  numerical 
answers  and  cause  these  predicted  answers  to  approach  or  deviate  from  the  experimental  results 
can  be  quickly  determined. 

Moreover,  once  these  numerical  solutions  are  available,  they  also  become  of  value  in  form¬ 
ing  a  foundation  for  subsequent,  more  complex  analytical  procedures.  This  approach  has  been 
adopted  in  the  numerical  analysis  of  the  cylindrical  shells.  In  Section  3  are  the  results  of  the  cal¬ 
culations  using  analytical  methods  based  on  linear-elastic,  small  displacement  theory.  Included 
are  comparisons  of  various  combinations  of  overpressure  and  positive  phase  durations. 

The  information  and  data  obtained  in  the  limited  experimental  portion  of  the  program  are 
presented  in  Appendices  I  and  II.  These  data  (involving  the  radial  displacements)  are  compared 
with  the  analytical  results  in  Section  4. 


Manuscript  released  by  the  authors  July  19^2  for  publication  as  an  ASD  Technical  Documentary 
Report . 
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2.  ANALYTICAL  PROCEDURES 


2. 1  Response  of  Cylindrical  Shells  to  Impulsive  Loadings 

2.1.1  Linear -Elastic,  Small- Deflection  Theory 

The  development  of  the  necessary  equations  for  the  analysis  of  cylindrical  shells  begins 
with  the  displacement  [Equation(I.  12)]  and  frequency  [Equation(I.  13)]  equations  in  Appendix  I.  We 
assume  that  the  cylinder's  end-closures  consist  of  elastic  media  which  offer  unequal,  elastic 
restraint  against  rotation.  For  generality,  we  define  these  elastic  restraints  as  Sj  at  x  =  0  and 
^2  ^  ^  where  S  is  the  stiffness  per  unit  length  of  the  medium  or  the  moment  required  to  rotate 

a  unit  length  of  the  medium  through  an  angle  of  one  radian. 

The  boundary  conditions  are  taken  as 


-S, 

I  ^  1  0  X 

X  =  0 


=  0 


M 


-  S- 


0W 


x<f>  “  ‘^2  ax 

X  =  L 


w  =0 

X  =  0 


w  =0 

X  =L 


X  =  0 


X  =L 


(2.1) 


The  first  two  equations  denote  the  condition  of  continuity  between  the  shell  and  the  end-closures. 

The  condition  that  the  radial  displacements  are  zero  at  the  shell's  boundaries  appears  reasonable, 
particularly  for  end-closures  (such  as  flat  plates)  which  offer  a  high  degree  of  transverse  restraint. 
Since 


Equations  (2.1)  may  be  rewritten  as  follows 


0^W  1 

0W 

0X^ 

-  Rl 

0X 

lx  =  0 

X  =  0 

^2 

0  W  1 

0W 

-  ^2 

0X 

lx  =  L 

X  =  L 

w| 

0 

:=  0 


|x  =  0 

w|  =0 

X  =L 


(2.2) 


2 


where 


Utilizing  Equations  (2.  2)  in  Equations  (I.  12)  (Appendix  I)  we  obtain  the  following  set  of 
linear,  homogeneous  equations 

2  •  0  -I 


2  =  » 

j 

y  x^c.  -  R,  y  x.c.  =  0 

y  j  j  1  Y  j  J 

Z  -  R2  Y,  =  0 

j  j 


(j  =  1,2,  3,4) 


Since  Xj  =  k,  =  -k,  X^  =  iK,  X^  =  -iK  (see  Appendix  I),  Equations  (2.3)  yield  the 
following  coefficient  determinant  which  must  vanish  in  order  for  these  equations  to  be  consistent. 


-(k  +  iRj) 


-(K  -  iRj) 


I  (K  -  R^)e'^  (K  +  ■  ^^'2^®'"*^  I 

Expanding  the  determinant,  we  have 

A  =  8i[2fc^  sinh  fc  sin  fc  +  k(R^  -  R2)(sinh  #c  cos  K.  ~  cosh  K  sin  k) 
+  RjR2(1  -  cosh  K  cos  fc)  =  0 


Given  the  appropriate  values  of  Rj  and  R^*  corresponding  values  of  k  may  be  determined, 

thereby  providing  the  necessary  values  for  =  k/L  for  the  frequency  equation 

[  Equation  (I.  13)]  . 


The  related  displacement  equations  are 

KofcC  (  zl  -  1 

’^mn^^’ ~  r  K  i7  cosh-j7  +  sinh  /c  cos  j 


’^mn^^’ -  TjKjL 


r  /  'tx 

KX  ^ 

i  /  k:x\1 

[rz  (sinh^ 

+  sin  — 

J 

)  +  r4(-cosh  -j;jJ 

cosh  K  cos 


-)  -  R1R2  |^r3(sinh-!^ 


(/cx 

CO"  — 


3 


/  1  \  ^2^  r  2  /  Kx  Kx  \ 

^  sin  k  sinh  —  +  sinh  K  sin  j 


+  kR 


1  h(' 


fCX  Kx 

cosh  -j-  -  cos 


)  +  r4  (sin-^^  -  sinh-!^jj 


+  ^kR^ 


cos  K  sinh 


Kx 


cosh  K  sinh 


KX 

cosh 


-  cos 


(2.7) 


'^mn(x>  <!>,  t) 


c 


cos  m(|)  sin 


r  .  •  u  •  u  • 

,  sin  K  sinh  -=—  +  sinh  K  sin  — 
^  L  L 


+  kRi 


r2  ^cosh  -  cos  +  r4  ^sin-!^  -  sinh  jj 


(‘ 


+  2kR2  ^cos  k  sinh -j- 


cosh  K  sin -  RjR^  j^r3  ^cosh 


-cos 


+  r. 


^  .  K.X 

Sin  *r —  -  sinh 

I  ^ 


(2.8) 


where 


Fj  -  2k  sin  K  +  KRj(sinh  K  -  cosh  K  +  sin  K  -  cos  k)  +  ZkR^  cos  k 
+  Rj  R2(sinh  k  -  cosh  k  +  sin  k  +  cos  k) 


^2  = 

sin  K  +  sinh  K 

II 

cosh  K  -  cos  K 

II 

cosh  K  +  cos  K 

Having  defined  the  relation  for  the  deflections,  it  then  beconaes  a  simple  matter  to  obtain  the 
accelerations  as  functions  of  time. 

For  the  case  where  the  elastic  restraints  at  either  end  of  the  cylinder  are  equal 
(Rl  =  R2  =  R),  Equation  (2.  5)  reduces  to 


A 


8i[K^  sinh  K  sin  K  +  R^(l 


cosh  K  cos  k)]  =  0 


(2.9) 


and  Equation  2,  8  becomes: 


4 


<t>.  t)  =  p-  cos  m4>  sin  to^nt 


^si 


in  K  sinh  -jj  +  sinh 


KX\ 

fc  sin^j 


-  (CR  j^r2  (cos^  -  cosh^j  +  r3  ^sinh-^^+  sin-^^^j 

-  ^r3  (cosh-!^  -  cos-!^)  +  ^sin— -  sinh  — )]} 


where 


=  Z/c  sin  fc  +  R(k  +  R)[sinh  n  -  cosh  K  +  sin  K  +  cos  K.] 

Similar  expressions  for  u  (x,  ct>.  t)  and  v^Jx,  <j),  t)  are  found  by  using  the  relations 


u^„{x,  (I),  t) 


K3  9wj^j^(x,  <|>,  t) 


^2  9Wjnn(x.<|>,t) 

tiK,  9<t) 


{2.  10) 


(2.11) 


If  the  cylinder's  end-closures  are  extremely  rigid  such  that  the  boundary  conditions 
may  be  assumed  fixed,  then  R  — ►  00  and  Equations  (2.  5)  and  (2.  8)  reduce  to 

A  =  1  -  cosh  K  cos  K  -  0 


(2.12) 


where 


Wmn(x.  tlJ.t)  =  i^T  cos  m4.  sinwj^jjt  j^r3  (cos -j;- 

(,  /cx 

sintr^  -  sin  —  jj 


r”  =  sinh  K  -  cosh  /c  +  sin  k  +  cos  /c 

1 


(2.13) 


For  the  condition  where  the  cylinder* s  end-closures  are  such  that  the  boundaries  are 
simply  supported,  then  R  =  0  and  Equations  (2.  5)  and  (2.8)  become 


A  =  sinh  K  sin  k  =  0 


'"'mn  ‘t’-  = 


C4  cos  m<|>  sin  [Ref-  Eq.  (2.  3)] 


(2.14) 


(2.15) 
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Equations  (2.  12)  and  (2.  14)  provide  the  extreme  values  of  K.  That  is,  for  the  cylindri¬ 
cal  shell  with  both  edges  fixed,  Equation  (2.  12)  gives 


K  =  Y  +  1]  TT  =  n*ir 


(n  =  1, 2, 3, ... ) 


(2.16) 


whereas  for  the  simply  supported  boundaries.  Equation  (2.  14)  gives 

K  =  mr  (n  =  1, 2,  3,  .  .  .  )  (2.  17) 

Similarly,  from  Equation  (2.9),  we  obtain  the  variation  of  /c  with  the  term  representing  the  end- 
closures*  stiffnesses.  This  is  shown  in  Figure  2.1. 

2.1.2  Linear- Elastic,  Large  Deflection  Theory^ 


In  order  to  avoid  the  prohibitive  difficulties  associated  with  a  detailed  analysis  and 
thereby  obtain  a  solution  that  will  lead  to  useful  information,  we  assume  that: 

(a)  The  elements  normal  to  the  middle  surface  remain  normal  and  unstretched 
(extensional  or  membrane  vibrations) 

(b)  The  cylinder  wall  is  thin  ^say  —  < 

(c)  The  longitudinal  (u)  and  tangential  (v)  displacement  components  are  small 
compared  to  the  radial  displacement  (w). 

For  the  equations  of  motion,  we  have 


SN 


aN 


xy 


ay 


=  phu 


aN 


xy 

ax 


phv 


(2.18) 


Following  the  format  of  Fung  and  Sechler^^*  we  define  the  strain  components  at  the  shell's 
middle  surface  as 


au 

i/ 

awr 

ax  ^ 

2  \ 

ax  / 

av 

aw\^ 

w 

ay 

2  \ 

.  ay  / 

a 

au 

a 

V  aw 

aw 

ay 

X  ^  ax 

ay 

(2.19) 


By  neglecting  the  inertia  terms  in  the  x  and  y  directions  (this  corresponds  to  a  first  order  approxi¬ 
mation  in  the  perturbation  procedure)  and  introducing  the  stress  function  F(x,  y)  where  Equa¬ 
tions  (2.  17)  are  satisfied  by 


a^F 


ay 


2  ~  ^xx’ 


a^F 

ax^ 


^Notation  in  this  section  same  as 


^yy’  ■  axay  " 


that  given  in  Appendix  I  except  where  noted. 
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.S/M/=>LE.  SUPPOaT 


FIGURK  2.  1.  VARIATION  OF  FREQUENCY  PARAMETER 
WITH  STIFFNESS  PARAMETER 
(See  Eq.2.  9) 


we  obtain  the  compatibility  equation 


Eh 


F/  9^w  \ 

^2 

0  W 

,,2 

0  W 

1 

9"w] 

L\.9x9y/  - 

9x2 

9y2 

“  a 

9x2  J 

and  the  equation  of  motion  for  the  radial  displacement 


(2.  20) 


4 

V  w 


1  w 


9^F 

9y2 


-  2 


„2 

0  W 


0^F 


0x0y  0x0y 


,2 

0  W 


0^F 


1 

+ 

a 


9^F 


0x' 


Y  +  p(x,  <|),  t) 


-  pwj 


(2.21) 


The  problem  now  becomes  one  of  assuming  a  form  for  the  radial  displacement,  using 
it  to  determine  the  stress  function  from  Equation  (2.  20)  which,  in  turn,  is  applied  to  Equa¬ 
tion  (2.  21).  The  latter  equation  suggests  that  the  radial  displacement  may  be  of  the  form 


w(x,  y,t)  =  f(t)  *  G(x,  y) 


where  G(x,  y)  denotes  the  mode  shape. 


The  most  appropriate  expression  is  one  which  is  sufficiently  complete  so  as  to  embrace 
all  the  significant  parameters  but  not  unduly  complex  so  as  to  complicate  the  analysis  and  intro¬ 
duce  refinements  that  are  not  commensurate  with  the  required  accuracy  of  the  final  numerical 
answers. 


As  an  example  of  a  simple  mode  shape,  we  take 

w(x,  y,  t)  =  f(t)[  sin  ax  cos  Py  +  sin  ax  cos  yy]  (2.  22) 

where  a  =  'rr/L,  p  =  1/a,  y  -  m/a,  and  m  denotes  the  number  of  circumferential  half  waves. 

Using  Equation  (2.  22)  in  Equation  (2.  20)  and  after  extensive  algebraic  manipulation,  we  arrive  at 
the  following  expression  for  the  stress  function 


where 


Eh[  f(ta^a^  , 

F  =  — ; - 7 -  [Kj  cos  2py  +  K2  cos  2yy  +  K3  cos  (p  -  7)y  +  K4  cos  (P  +  7)y 


+  cos  2ax  +  cos  (p  -  7)/  cos  2ax  +  cos  (p  +  7)y  cos  2ax 

+  Kg  sin  ax  cos  py  +  sin  ax  cos  yy] 


(2.23) 


K, 


8p^ 


- 5* 

4  (p  +  7)2 


K 


4'(P  -  y) 


”  [4a2  +  (P  +  7)2]  2 


K 


2  -  -  2 
87 


Kp 


8a4 


Ko  = 


f(t)a(a^  +  p2)2 


K3  =  - 


(P  -  7) 


Kp  = 


4'  (P  +  g)*^ 


[4a^  +  (p  - 


Ko  = 


4^ 


f(t)a{a^  +  yY 
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For  a  radial  displacement  equation  of  the  form 

w(x,  y,  t)  =  f(t)[sin  ax  cos  |3y  +  sin  ax  cos  yy  +  ^  sin  ax] 
the  stress  function  is  as  follows 


F  = 


Eh[f(t)] 


[Kj  cos  2Py  +  cos  2yy  +  K3  cos  (p  -  y)y  +  cos  (P  +  Y)y 


+  cos  Zyx  +  cos  (j3  -  y)y  cos  Zax  +  cos  (P  +  7)y  cos  2ax 
+  Kg  cos  2ax  cos  Py  +  cos  Zax  cos  yy  +  cos  py  +  cos  yy 

+  sin  ax  cos  py  +  Kjg  sin  ax  cos  yy  +  sin  ax] 

where  Kj  through  Ky,  are  the  same  as  those  given  for  Equation  (Z.  Z3)  (except 

that  ^  =  1)  and 


K 


Z^p^ 


^  ^  (4a^  +  P^^ 


K 


Z£ 

10  -  „Z 


K 


Z^y 


^  ~  (4a^  +  P^)^ 


K 


Z$ 
11  =  Z 


As  a  final  example,  for 

w(x,  4>,  t)  =  f(t)[A  sin  ax  cos  py  +  sin  ax  cos  7y  +  X  sin^  ax] 
the  stress  function  becomes 

Z  Z 

p  =  cos  zpy  +  cos  Zyy  +  Kg  cos  (P  -  y)y  +  K^  cos  (P  +  y)y 

+  Kg  cos  Zax  +  Kg  cos  (p  -  y)y  cos  Zax  +  K^  cos  (P  +  y)y  cos  Zax 
+  Kg  sin  ax  cos  Py  +  K^^  sin  ax  cos  yy  +  K^^  sin  3ax  cos  Py 
+  K^j  sin  3ax  cos  yy] 

where  Kj  through  K^,  Kg,  and  K^  are  identical  to  those  for  -Equation  (Z.  Z3)  (for  =  1)  and 


y  +  (m  - 


z  4x 


Kg  = 


af(t) 


8a^ 


^  [1  ■ay^x^(l)]4 
f(t)a(a^  + 


(2.  2-4) 


(2.25) 


(2.  26) 


(2.27) 
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_  [1  -  ap\f(t)]4q/ 

^  '  f(t)a(a2  +  p2)2 


K 


XT 


10  ■  2  ^  2.2 

(9a  +  y  ) 


K 


■  (9a2  +  p¥ 


The  Galerkin  procedure  may  now  be  used  to  find  the  differential  equation  for  the  time 
function  f(t),  thereby  providing  not  only  the  deflection,  but  also  the  acceleration  as  a  function  of 
time.  The  Galerkin  equation  will  have  the  form 


a^F 


axBy 


+ 


a^w  a^F  a^F 

0y2  9x^  ^  9x^ 


(2.  28) 


where,  for  example,  the  stress  functions  (F)  derived  from  the  radial  displacement  expression 
w(x,  <t),  t)  would  be  those  from  Equations  (2.  22)  and  (2.  23),  (2,  24)  and  (2.25),  or  (2.  26)  and  (2.  27). 

In  order  to  obtain  a  solution  for  Equation  (2.  28),  it  is  first  necessary  to  specify  the 
external  loading  conditions,  p(x,  {j),t).  One  satisfactory  relation  is  the  following 


p(x,  4),  t) 


( 


Po  +  Pi 


(2.29) 


where,  for  example, 

T(t)  =  e-^t 


or 


T(t) 


ct 


The  differential  equation  for  the  time  function  is  of  the  form 
J  +  +  AjJ  =  A4T(t) 


(2.30) 


where  J  =  f(t)/a  and  are  coefficients  reflecting  the  particular  characteristics  of  the 

mode  shape  and  load  function.  As  a  particular  example,  we  take  L  =  36  inches,  A  =  6  inches, 
h  =  .  036  inch,  E  =  2.  9  X  10^  psi,  v  =  0.  25  and  p  =  24.  7  X  10"^  Ib-sec^/in^  as  the  properties  of 
the  cylindrical  shell.  For  the  mode  shape  defined  in  Equation  (2.  26)  (with  x  =  0*  or  0, 

and  n  =  2)  and  the  loading  defined  as 


p(x,  4),  t)  =  Pq(1  +  B  cos  Py)e 


-6t 


the  numerical  values  of  the  coefficients  are  as  follows: 
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For  A=  1: 


pp  =  70  psi,  B  =  4/7  and  6  =  2.77  X  lO'^ 

Ai  =  23.4,  A2  =  -14.3,  A3  =  5.8,  A4  =  3.72  X  10'^ 

For  A=  0:  Pp  =  70  psi,  B  =  4/7,  and  6  =  2.  77  X  10 

Aj  =  10.4,  A2  =  -1.1,  A3  =  0.3,  A4=  1.65X10’^ 

The  variation  of  J  as  a  function  of  T(t)  is  given  in  Figure  2.  2.  Accordingly,  the  maximum  radial 
displacements  for  A=  1  and  A  =  0  are  respectively 


w 


max 


0.0Z2" 


w 

max 


0.  11" 
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7~  C^)(^s<&c,) 

FIGURE  2.  2.  VARIATION  OF  TIME  FUNCTION  WITH  T(t) 
{See  Eq.  2.  30) 


12 


2.  2  Dynamic  Response  of  Circular  Plates 
2.  2.  1  Large  Amplitudes 

LIST  OF  SYMBOLS 


a 

D 

E 

f,  F,  G 
h 

•^n>  ^n»  ^n>  ^n 
N 

m 

Nr,  Nt 

P 


r,  0 


r 


t 

I 


V 


w 


radius  of  circular  plate  (in.  ) 

- - ,  flexural  rigidity  of  plate  (lb  in. ) 

12(1  -  v2) 

modulus  of  elasticity  of  plate  (Ib/in^) 
functions  of  time 
thickness  of  plate  (in.  ) 

Bessel  functions  of  order  n 
(Nr  +  Nt)/(1  +  v)  (lb  in"^) 
buckling  load  of  circular  plate  (lb  in"^) 

stress  resultants  in  the  radial  and  tangential  directions,  respectively  (lb  in"^) 
natural  frequency  of  plate  (sec"^) 
cylindrical  coordinates 

£ 

a 

time  (sec) 


th 

a2 


2(1  -  v^)p* 


velocity  (in.  sec“^) 


a^  / 12(1  -  v^)p* 

''h2  V  E 


lateral  deflection  of  plate  (in.  ) 


w 

h 


mass  per  unit  area  of  plate  (lb  sec^  in 
mass  density  of  plate  material  (lb  sec 


=  (V^)^  =  biharmonic  operator 


2 in~4) 


+ 


I _ 

r2  9  02/ 


2 


It  is  well  known  that  when  the  lateral  deflections  are  large  the  behavior  of  plates  is  gov¬ 
erned  by  two  coupled  nonlinear  partial  differential  equations  usually  credited  to  von  K^rm^n. 

The  dynamic  analogue  of  these  equations  has  been  derived  by  G.  Hermann,  Chu  and 

Hermann(2* have  calculated  the  fundamental  frequencies  of  rectangular  plates  at  large  amplitudes 
by  an  energy  procedure. 
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However,  a  theoretically  exact  solution  of  the  dynamical  equations  even  for  the  case  of 
free  vibrations  is  nonexistent.  The  case  of  forced  vibrations  had  apparently  not  been  studied  at  all. 

It  is  to  be  realized  that,  owing  to  the  nonlinearity  of  the  equations,  the  free  and  forced  vibration 
problems  cannot  be  solved  separately  and  superposed,  neither  can  advantage  be  taken  of  the  concept 
of  normal  modes  as  is  usual  in  the  linear  case. 

In  1955,  Berger(^*  proposed  a  set  of  equations  which  may  be  termed  a  simplified  version 
of  the  von  K^rm^n  system.  The  merit  of  his  method  lies  in  decoupling  the  two  equations,  so  that  one 
of  them  assumes  a  quasi- linear  form  and  can  readily  be  integrated.  He  showed  further  that  his  equa¬ 
tions  yielded  solutions  remarkably  close  to  those  obtained  by  more  elaborate  procedures  using  von 
K^rm^n*s  equations,  at  least  as  regards  deflections.  Berger,  it  must  be  remarked,  confined  him¬ 
self  to  the  static  case, 

Nash  and  Modeer^^*  extended  the  Berger  equations  to  the  dynamic  case  for  rectangular 
plates  and  showed  that  the  free  vibration  problem  can  be  treated  to  yield  results  which  are  reason¬ 
ably  close  to  the  results  obtained  by  Chu  and  Hermann.  (2- 4) 

In  what  follows,  the  dynamical  equations  a  la  Berger  are  extended  to  circular  plates.  By 
using  an  approximate  step-by-step  procedure, it  is  shown  that  not  only  the  free  vibration  but  also  the 
forced  vibration  problem  may  be  solved  using  the  elementary  concept  of  normal  modes.  It  is  neces¬ 
sary  to  add  that  although  the  basic  Berger  formulation  has  a  quasi- linear  form  for  the  static  case,  it 
is  nonlinear  for  the  dynamic  case  and  cannot  be  solved  with  comparable  ease. 


Z.2.1.1  Equations  of  Motion 

For  the  dynamical  case,  the  Berger  equations  for  circular  plates  may  be  written  in  the 

form 

v^w  -  a2v2w  +£- 

D  D 

g^h^  =  +  i.  f  9w \  ^  H  +  i.  +  1  /  9w\  ^ 

IZ  3r  Z  \  9r  /  r  r  90  ^r^  \90/ 

in  which  =  (V^)^  is  the  biharmonic  operator,  p  is  the  mass  of  the  plate  per  unit  of  area,  w  is  the 

Eh^ 

lateral  deflection,  u  and  v  the  radial  and  tangential  displacements,  D  =  -  is  the  flexural 

1Z(1  -  vZ) 

rigidity  of  the  plate,  h  its  thickness,  E  is  the  elastic  modiilus,  v  is  Poisson’s  ratio  and  f  a  given 
forcing  function. 


(2.31) 


The  radial  and  tangential  strains  and  may  be  expressed  in  terms  of  the  stress  result¬ 
ants  Nj.  and  as  follows 


er  =  (Nr  -  vNt) 

et  =  (Nt  -  vNj.) 


(2.32) 


It  follows  from  (Z.  3Z)  and  the  second  of  (Z.  31)  that 


2u2 

ah 

IZ 


1  -  V 


Nh^ 


(2.33) 


In  the  Berger  formulation,  a  is  constant  throughout  the  plate  and  thus  N,  which  has  the 
same  dimension  as  a  stress  resultant,  is  also  constant  throughout  the  plate.  [The  simplicity  of  the 
Berger  equations  apparently  derives  from  the  basic  assumption  that  (Nj.  +  N^)  is  constant.  ] 
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From  Equations  (2.  33)  and  the  second  of  (2.  31)  one  finds 

+  -L  /  \  ^  i  ^  ^  \ 

12D  3r  2  \  9r  /  r  r  BO  2r^  \99/ 

In  view  of  the  assumed  constancy  of  N  with  respect  to  r  and  0,  one  may  multiply  by  r  dr  d0  and  inte¬ 
grate  over  the  plate  to  find 


If  the  boundary  conditions  are  such  that  u  and  v  vanish  at  the  boundaries,  the  first  two  integrals  on  the 
right-hand  side  are  evanescent  and  there  must  follow 


Nh^ira^ 

12D 


2 


2tt  a 

II 


0  0 


rdr  de 


using  (2.  33)  the  first  of  (2.  31)  may  be  written 


V'^ 


N  v2w  +  £-  —  = 

D  D  at2  D 


(2.34) 


The  problem  now  reduces  to  the  integration  of  the  system  (2.34). 


2.  2. 1.2  Integration  Procedure 

It  is  seen  from  the  first  of  (2.  34)  that  N  is  a  function  of  time  alone  and  that  the  second  of 
(2.34)  is  effectively  nonlinear.  The  problem  of  finding  the  dynamic  response  of  such  a  system  cannot 
be  solved  by  superposition  of  the  separate  solutions  of  free  vibrations  and  of  the  steady  state  prob¬ 
lem. 


A  step-by-step  integration  procedure  is  here  proposed  for  the  integration  of  the  system 
(2.  34).  The  assumption  is  mac^  that  N,  which  is  proportional  to  the  sum  of  the  membrane  tensions 
(Nj.  +  N^),  remains  constant  at  N  for  a  sufficiently  short  interval  of  time ^  .  It  is  evident  that  the 
second  of  (2.  34)  assumes  a  linear  form  and  may  be  readily  integrated.  N  is  evaluated  from  the  first 
of  (2.  34)  at  the  beginning  of  each  time  interval.  The  nonlinear  problem  is  thus  reduced  to  a  series  of 
lineal*  problems  in  time. 

Even  this  simplification,  it  may  be  noted^  still  involves  a  prodigious  amount  of  labor.  For, 
since  the  normal  modes  of  vibration  are  a  function  N,  they  are  a  function  of  time  also.  Furthermore, 
the  forcing  function  f(r,  0,  t)  has  to  be  expanded  in  a  series  of  the  normal  modes  appropriate  to  that 
time  interval. 

Thus,  the  general  procedure  is  as  follows:  Starting  with  the  given  initial  conditions,  eval¬ 
uate  N  from  the  first  of  (2. 34).  Use  this  value  of  N  to  solve  the  homogeneous  part  of  the  second  of 
(2.34)  and  determine  the  natural  frequencies.  Expand  the  forcing  function  in  terms  of  the  normal 
modes  and  determine  a  particular  solution  of  the  second  of  (2.  34).  Use  the  initial  conditions  at  time 
zero  to  determine  the  deflection  and  velocity  at  the  end  of  the  time  interval  A.  Use  the  latter  as  the 
initial  condition  for  the  second  time  interval  A.  Repeat  this  procedure  up  to  any  desired  time  t  =  sA. 

(Some  of  the  tediousness  of  this  procedure  may  be  reduced  if  one  ignores  the  change  in  the 
normal  functions  at  each  time  interval,  using  throughout  instead  the  functions  appropriate  to  the  initial 
interval,  but  taking  account  of  the  change  in  the  natural  frequencies  at  each  interval.  ) 
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Returning  now  to  the  second  of  (2.  34),  let 


w  =  R(r)e(0)  sinpt  (2.  35) 

where  R  is  a  function  of  r  alone  and  0  is  a  function  of  6  alone.  On  substitution  into  the  homogeneous 
part  of  the  second  of  (2.  34)  and  dividing  out  the  common  factor  sinpt  there  follows 

v^(Re)  -  ^  v^(Re)  -  ^  Re  =  0 

which  may  be  written  in  the  form 

, 2 \  /  ^  q2 


(2.36) 


where 


.2  _ 


Na^ 

[(■ 

1/2 

2D 

n2  / 

Na2 

[(■ 

+  iE!£2\ 

1/2 

2D 

) 

a2  = 

Na^ 

D 

■] 

3 


(2.37) 


=  /E  ^ 

V  P 


(2.38) 

(2.39) 


Equation  (2.  36)  shows  that  the  complete  solution  may  be  obtained  by  adding  together,  with 
appropriate  arbitrary  constants,  the  solutions  of  the  two  equations 


and 


Re  =  0 


RG  =  0 


(2.40) 


On  letting 
*  0  =  cos(n0  -  7j^) 

Equations  (2.40)  yield  two  ordinary  differential  equations  for  the  determination  of  R 

isa,  (si.  R.o 

dr2  r  dr  V  a2  / 

iaa.  r.« 

dr2  r  dr  V  a2  r2  / 

The  solutions  of  (2.42)  are,  respectively 

(f) 


(2.41) 


R  =  AJ, 


(2.42) 


and 
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in  which  and  denote  Bessel  functions  of  the  order  n  of  the  first  and  second  kinds,  respectively. 
Ill  and  Kn  are  the  modified  Bessel  functions  of  the  first  and  second  kinds,  and  A,  B,  C,  and  D  are 
arbitrary  constants.  (The  notation  for  the  Bessel  functions  is  that  used  in  Reference  2.  7. ) 

In  the  case  of  a  solid  plate,  Yn  and  Kn  are  inadmissible  because  of  the  singularity  at  r  =  0, 
and  the  general  solution  of  the  system  (2.42)  is  thus 


R  -  AnJn  +  Bnin  ( 

The  general  solution  of  the  homogeneous  part  of  the  second  of  (2.  34)  may  now  be  written 
w=|^AnJn^^)  |cosn0  +  CnSinnej.  I^sinpnt  +  Dncospntj.  (2.44) 

The  frequency  equation  and  the  ratio  B^/ may  be  obtained  from  the  boundary  conditions. 


In  the  case  of  a  simply  supported  plate,  one  has 


w  =  0 

9  ^w  f  ^  9w  9  \ 

Br  ^ 

and  in  the  case  of  a  clamped  plate 


(2.45) 


§7  =  " 


at  r 


(2.46) 


Substitution  of  (2.44)  into  (2.45)  yields 

Bn  _  Jn(a) 

^  ■  ■  In(P) 

•^n  +  1^“^  o  +  l(P)  _ 

Jn(<i)  In(P)  i  - 


(2.47) 


and  the  substitution  of  (2.44)  into  (2.46)  leads  to 


^n  +  1^^^^ 
Jn('i) 


+  P 


InH-l(P) 

MP) 


=  0 


> 


J 


(2.48) 


The  second  of  (2.47)  is  the  frequency  equation  for  the  simply  supported  plate  and  the  secofid  of  (2. 48) 
for  the  clamped  plate.  These  equations,  together  with  (2.38),  are  sufficient  for  the  determination  of 
the  natural  frequencies  p^  for  given  values  of  N. 


17 


Defining 


N*  ^ 


4.  2D 

—  with  V  =  0.3,  for  a  simply  supported  plate 

a2 

14.68D 


for  a  clamped  plate 


(2.49) 


the  frequency  equation  for  the  simply  supported  plate  may  be  written 


Jn  +  l(-)  .  „  In  +  l(P)  p2+a2 

a — r-T— ; —  +  6 — =— rr: - =  - 


Jn(“) 

=  4.24) 


V  =  0.3 


UP) 


1  -  V 


(2.50) 


and  that  for  the  clamped  plate 

■^n  +  1^*^)  ,  o  ^  _  « 

“  Jn(a)  P  In(P) 

p2  -  =  14.684) 


(2.51) 


Equations  (2.  50)  and  (2.  51)  may  be  solved  for  various  values  of  the  parameter  ^  and  the 
values  of  a  and  p  so  derived  may  be  tabulated  or  graphed  as  functions  of  <[). 


It  will  be  noted  from  (2.  45)  and  (2.  46)  that  the  ratio 


fnj  ^  Jn(qj) 
■^nj  In(Pj) 


whether  the  plate  is 


simply  supported  or  clamped. 

Letting 

one  may  write  the  deflection  function  (2.44),  corresponding  to  the  (nj)^^  frequency  as  follows 

=  {A^jCOsnO  +  Bn^jSinnO  }  (sinpj^jt  +  D^j  cos  p^jt  |  (2.53) 


"nj 


in  which  A^j,  and  D^j  are  new  arbitrary  constants,  to  be  determined  from  the  initial  conditions. 
It  may  be  seen  that  (2.  53)  may  be  written  in  the  form 


''nj 


^nj  '{^nj(t)  cos  n0  +  Bnj(t)  sinns]- 


where 


Anj(t)  =  Anj  -^sinpnjt  +  Dnj  cospnjt^ 
Bnj(t)  =  Bnj  •(sinpnjt  +  Dnj  cospnjtj- 


(2.54) 


(2.55) 


Any  arbitrary  deflection  satisfying  the  boundary  conditions,  and  therefore  the  solution  of  the  homo¬ 
geneous  part  of  the  first  of  (2.  34)  appears  finally  in  the  form 
OO  00  00  oo 

w^  =  S  S  RnjAnj(t)  cos  n0  +  E  E  RnjBnj(t)  sinn  (2.56) 

n=0j=l  n=:lj=l 

where  the  superscript  "c”  on  the  left  denotes  "complementary  function.  " 
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2.2.  1.3  Forced  Vibration 


In  seeking  a  particular  solution  of  the  first  of  (2.  34),  it  will  be  supposed  that  the  forcing 
function  f(r,  0,t)  can  be  expanded  in  terms  of  the  principal  modes  in  the  form 

f(r,  0,  t)  =  F(t)  SSR„;(rj^^  cos  n0  +  s^i  sin  n0)  (2.  57) 

n  j  ^ 

A  particular  solution  of  (2.  34)  is  sought  in  the  form 

wP  =  SSRjij  ’l^anj(t)  cos  n0  +  bn.j(t)  sinn0  (2.  58) 

in  which  a^j  and  b^j  are  undetermined  time  functions. 

Inasmuch  as  R^j  cosn0  and  Rj^j  sinn0  satisfy  the  homogeneous  part  of  the  first  of  (2.  34) 
there  must  follow 


(cos  n0  \  /  cos  n 

sinn0  /  \  sinn( 


Z 

PnjP 


COS  n0 ' 


sinn0 


(2.59) 


Substituting  (2.  58)  into  the  first  of  (2.  34)  and  making  use  of  (2.  59)  one  finds 


anj  +  Pnjanj 


F(t)f, 


Hi 


bnj  +  Pnjt>nj  “  ' 


P 

F(t)s 


iLl 


(2.60) 


} 


Let  the  particular  solutions  of  these  equations  be 
anj  =  injG(t) 

=  hnjG{t) 

From  (2.  58)  and  (2.  61),  one  gets 

wP  =  SSRjjj  (gjjj  cos  nO  +  hj^j  sinn0)G(t) 

The  complete  solution  of  the  first  of  (2.  34)  may  now  be  obtained  by  adding  (2.  56)  and  (2.  62) 

oo  oo 

w  =  w~  +  w''  =  z;  i  K„;lA„;lt)  +  g„;Utt)l  cosn0  +  S  S  R„ ; [ B„ ^ (t)  +  h„ iG(t)]  sin nd 


(2.61) 


(2.62) 


oo  oo 
.c  +  wP  =  S 
n  = 


0  .Sj  R^j[A^j(t)  +  g,jG(t)] 


S  S  Rnj[B„j(t)+h^jG(t)] 
n=lj=l 

(2.63) 


Let  the.  initial  conditions  be 
w  =  <|)(l 


w  =  <|)(r)4i(0) 


at  t 


These  can  always  be  expanded  in  Fourier- Bessel  series  of  the  form 


CO  00 
w  =  S  2 
j=l n=0 

00  00 


oo  oo 

RnjHnjcosn0  +2  2  RnjLnjsinn0 

j= 1 n= 1 


oo  oo 


2  2  RnjMnjcosn0+  2  2  RnjNnj  sinn0 

j=ln=0  j=ln=0 


(2.64) 
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where 


2u  a 

XHnj  =  J  J  <t>(r)4/(0)Rnj  cos  n0rdrd0 
0  0 


2u  a 

XLnj  =  j  J  <t)(r)4»(0)Rnj  sinner  dr  d0 
0  0 

2ir  a 

XMnj  =  J  J  cos  n0  rdrdO 

0  0 


n  ^  0 


XNnj  =  //  x(r)|jL(6)Rnj  sinnOrdrdO  n  ^  0 


(2.65) 


in  which 


I 


\  =  Ztt  I  Rnj^ 


Rnjr  dr 


n  =  0 


n  ^  0 


From  (2.63)  and  (2.64),  using  (2.65),  one  finds 

OO  CD  OO  CD 

w=  S  S  RniUniCOsn0+  S  S  RniW^iSinnO 

n=0  j=l  n=l  j=l 

CD  CD  OO  CD 

w=  2  2  RniVni  cos  n0  +  2  2  RniZnisinn© 

n=  0  j=  1  n=  1  j=  1 


in  which 


(2. 66) 


^  {*^nj  -  gnjG(O)}  cosp^.t  +  g^jG(t) 

’^nj  =  {^nj  -  Pnj*  "  Pnj  {*^nj  '  +  gnjG(t) 

^nj  +  {^nj  "  }  cos  +  h^jG(t) 

=  {^nj  -  hnjG(O)}  cosp^^t  -  p^^  [l^.  -  h^jG(O)  }  sinp^^t  +  h^jG(t) 

The  deflection  is  now  completely  known  and  the  acceleration  may  also  be  determined.  It 
is  to  be  remembered,  however,  that  this  solution  is  valid  only  for  the  time  integral  (t  -  r)  during 
which  the  parameter  N  is  supposed  to  remain  constant.  In  Equation  (2.66),  therefore,  the  time  vari¬ 
able  must  be  (t  -  t)  instead  of  t.  The  parameter  N  is  determined  by  the  first  of  Equations  (2.  34)  with 
t  taken  at  the  beginning  of  the  time  interval  (t  -  r).  Thus 

Ztt  a 

N=-^  r  r 

Tra2h2  J  J 

0  0 


rdrde 


(2.67) 


t  =  T 
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Using  (2.  63)  and  (2.  66)  and  (2.  50)  and  2.  51),  Equation  (2.  67)  may  be  written 

RJ>jR{)kHojHok}‘rdr  + 


r2  „2  _  6 

^  h2 


[a. 

^  r  oo  GO 

0 

/  (  2  Z  S  RiijRlik(HnjHnk  +  I^jLnk) j  rdr  + 

Qln=lj=lk=l  J 

a 

^  r  CO  C»  CO  1  dr 

/  4  S  S  S  n^RnjRnkCHnjHnk  +  Imj^nk)  ^  — 

•/Ln=lj=lk=l  J  ^ 


(2.68) 


in  which  Ri^  =  — ; — - 
J  dr 

The  integrals  in  (2.68)  cannot,  unfortunately,  be  readily  obtained  in  closed  form,  and 
recourse  must  be  had  to  numerical  integration. 

2.  2. 1.4  Summary  of  Integration  Procedure 

In  numerical  integration  it  is  advantageous  to  use  nondimens ional  quantities.  These  will  be 
distinguished  by  placing  a  bar  over  the  corresponding  symbols.  Let 

—  w 
■  h 


1=  ^ 

a2  V  12(1  -  v2)p* 

where  p*  is  the  mass  density  of  material  of  the  plate 


_  r 
r  =  — 
a 


Let  s  denote  the  number  of  intervals  A  from  time  zero.  Then  the  deflections  and  velocities  at  the  end 
of  the  s^^  interval  are  given  by 


Z  Z  Rjii  Uni  cosn0  +  Z  Z  R^i  W^i  sinn6 
j=ln=0  j=ln=l 


00  oo 


oo  oo 


We  —  Z  Z  Rni  ^ni  cos  n0  +  Z  Z  Rni  ^ni  sinn0 
j=l  n=0  j=l  n=l 


>  s  >  1 


(2.69) 


The  Rnj  are  defined  by  Equation  (2.  52)  and: 
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Js 

+  {^nj^  -  gnj^G[A(s  -  1)]}  cos  (Pnjg^)  +  gnjg^t^®! 

'  {“nj.  -  -  ■"}  ' 

-  I’nj,  {»aj,  -  S„j^G[A(s  -  1)|}  .in  (p„j^A)  +  g„j^Gl  A=1 

*«i.  "fir  {””i.  ■  •  "’}  •“"’Aj/"  * 

“Js 

+  {^njg  “ 

Znjg  =  {Nnjg  -  hnjgGlA(s  -  1)]}  cos  (p„j^A)  - 

-  Pnjg  {^nj  -  hnj  G[A(s  -  1)]}  sin(p„j  A)  +  G[As] 


s>l  (2.70) 


cx> 


Hni  =  S 

k  = 


2  j  /  Unkg  .  1  Rnkg  .  1 


Rnjgrdr  /  J  R^j^rdr,  n  >  0,  s>2 


oo 


4  =  S  r  Vnk  ,  ^nk 

k  =  1  n  1  “  1 


Rnj^  rdr  j  j  ^njg  rdr,  n  >  0,  s  >  2 


L-njg  "  J  ^  ^*^js  -  1  -  1  j /  ^njg  rdr,  n  >  0,  s  >  2 

Njij^  =  ^  ^  ^nkg  _  J  _  1  ^jg  rdr  j ^  rdr,  n  >  0,  s  >  2 

and  Hnj^,  defined  by  the  initial  conditions. 


(2.71) 


The  natural  frequencies  are  given  by 


p  .  =  6  .  a  •  s  >  1 

Pnjg  Pnjg  njg 


(2.72) 
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's  s 


/•¥ 

00 

00 

=  12/1 

s 

s 

^oj 

0  ^ 

j=l k=l 

*  00 

00 

00 

+  6  /  i 

S 

s 

S 

0  ^ 

.n=l j=l 

k=  1 

r  °° 

00 

00 

+  6  /  \ 

s 

2 

0  ' 

Ln=l j=l 

k=  1 

>  S  >  1 


(2.73) 


with 


°Js  dF~ 

Some  of  the  above  relations  may  be  simplified  by  noting  that  if  A  is  taken  sufficiently  small 
.  I  ^  ^nj 

.1 

/  Rnk-  .  1  Rnjg  »  0  j 

0 

Equations  (2.  71)  then  become 


‘'S 

-  1 

n 

> 

0 

Mni 

Vni 

n 

> 

0 

■^Js  -  1 

s  >  2 

(2.74) 

Lnjg 

tfi 

^-Js  -  1 

n 

> 

0 

Nnj^ 

-  1 

n 

> 

0 

Equation  (Z.  73)  simplifies  to 

Pi  -  .§  •  12  EJSojkHoj.  H„k,  +  <■  JSfInjkIHnj,  H„k,  +  L„k^)  + 


nj 


where 


Sojk  / 
0 

lujk"/ 


^jk 


R  .  R  ,  4r 
*^Jo  ""^o  7 


If  the  right-hand  side  of  (2.  75)  is  evaluated  for  any  value  of  s,  then  setting 
p^_ag  =  4.2cj>  for  simply  supported  plate 

=  14.  68  <j)  for  clamped  plate 


(2.75) 


(2.76) 
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one  may  determine  ()>  and  thereby  the  values  of  a^j  and  pnj  these  have  been  tabulated  for  various 
values  of  (f),  [See  note  under  Equation  (2.  51).  ]  Equation  (2.  72)  may  then  be  used  to  determine  the 
nondimens ional  frequencies. 

2.  2. 1.5  Numerical  Example 

The  response  of  a  circular  plate  to  a  loading  pulse  will  be  investigated  by  the  method 
described  on  the  preceding  page. 

A  blast  of  short  duration  is  best  treated  as  an  impulse.  Assuming  that  the  blast  is  deliv¬ 
ered  uniformly  over  the  area  of  the  plate,  let 


ira^  j*  f(t)  dt  =  I 

0  ' 

6  being  the  time  at  which  the  pulse  decays  to  zero. 

By  the  momentum  principle,  the  uniform  starting  velocity  is 


(2.77) 


where  M  is  the  total  mass  of  the  plate. 


The  plate  under  consideration  has  the  following  geometric  and  material  properties 
a  =  6".  E  =  30X10^psi,  p*  =  7.  34  X  lO""*  lb  sec^  in"^ 


Thus, 


h2 


0.  45  sec  in 


The  initial  velocity  Vq.  calculated  from  the  characteristics  of  the  charge  (used  in  an  experiment,  the 
details  of  which  are  omitted  here)  and  relation  (2.  77)  above  gives  an  estimated  initial  velocity  of 
893  in.  sec"^.  Thus,  the  nondimens  ional  initial  velocity  is 


Vo  =  893  X  0.45  «  400 


The  response  of  the  plate  will  be  analyzed  on  the  assumption  that  the  plate  receives  a  uniform  velocity 
of  400  while  at  rest. 


The  parameter 


a2 

Thus 

t  =  6l . 6  t 

This  gives  the  relation  between  nondimensional  time  t  and  true  time  t.  The  nondimensional  time 
interval  A  will  be  taken  as  0.  06l6.  This  corresponds  to  a  true  time  interval  of  one  millisecond. 

The  basic  equations  given  in  Section  2.  2.  1.4  can  be  simplified  considerably  because  of 
axial  symmetry.  In  what  follows,  the  bars  over  the  nondimensional  quantities  will  be  dropped.  No 
confusion  can  arise  as  all  quantities  henceforth  are  nondimensional. 


12(1  -  v2)p* 


61. 6 
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Equations  (2.  69)  become 


Ws  =  Rojj  Uojg 

Ws  =  Rojj  Vojg 


s  >  1 


From  Equations  (2.  70),  (2.  71)  and  (2.  74) 

rj  TT  -  1  .  , 

^oj  ^oi  ,  (p. 

Js  ^Js  -  1  P„j 


°Js  -  1 


M^.  =  V„.  =  M„. 

OJs  °Js  -  1  °Js  -  1 


OJf 


jA)  +  ^  coMPoj_  _  A) 


cos  ^  A)  -  p„j^  _  _  ^  sin  (p„j_  _  ^  A) 


The  initial  conditions  give 


Hojj  -  0 


'of 


^oji  "  ^ 


Rr,,-.  rdr  /j  rdr 


^oji 


where 


with 


Thus 


Jo(^j]^ ) 

Roj  =  r)  -  I^(^.  ^  Io(“i,  r) 


Jl' 


aj^  =  2,  222  for  j  =  1 


=  5, 452  for  j  =  2 

=«f. 


M 


U„;  = 


oj, 


°Jl  P, 


oji 


sin  (PojA) 


V^i  =  cos  (pojA) 


OJl  OJl 
From  Equation  (2.  72) 


p  .  =  a  .  0  . 

Pojs  ojg  Pojg 


s  >  1 


Using  j  =  1,  2  only.  Equation  (2.  73)  becomes 


Assuming  that  the  plate  is  simply  supported 

Thus,  <j)g  is  known,  and  and  pg  may  be  read  from  charts  already  prepared. 


(2.78) 

f  s  >  Z 

(2.79) 


(2.80) 

(2.81) 

(2.  82) 
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2.2.2  Plastic  Collapse 


LIST  OF  SYMBOLS 

V  initial  velocity  given  to  membrane,  a  constant 

U  radial  velocity  of  flat  portion  of  membrane,  positive  outward 

R  radius  of  flat  portion  of  membrane 

R  radial  velocity  of  bending  wave 

a  initial  radius  of  membrane 

h  initial  thickness  of  membrane 

p  mass  density  of  membrane  material 

t  time 

r  radius  of  position  of  any  particle  in  the  flat  portion  of  membrane  at  any  time 

Tq  initial  radius  corresponding  to  r 

H  thickness  of  membrane  at  any  time 

(Tj.  radial  stress 

O’ Q  circumferential  stress 

O'  yield  stress  in  tension 

Cj.,  natural  (logarithmic)  strains  in  the  radial  and  circumferential  directions 

distance  of  central  flat  portion  from  initial  plane  of  the  diaphragm 
Poisson*s  ratio 
modulus  of  elasticity 
=  R/a 
=  U/v 
H/h 
=  (r/p 

=  vt/a=  Z/a 
=  40-  /  pv^ 

=  [(4o-/pv2)cr(l  -  v)/E]  1/2 

In  investigating  this  problem  it  was  discovered  that  there  are  certain  existing  solutions  to 
corresponding  problems  in  underwater  shock  which  have  been  satisfactorily  verified  experimentally. 
While  the  problems  of  blast  effects  in  air  or  of  soft  radiation  are  by  no  means  the  same  as  the  effects 
of  blast  loads  underwater,  they  are  sufficiently  similar  to  warrant  the  supposition  that,  as  a  first 
approximation,  the  two  problems  may  be  studied  by  identical  procedures. 


Z 

V 

E 

P 

c^ 

t 

6 

K 
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The  solution  to  the  plastic  collapse  of  a  diaphragm  proposed  by  HudsonC^*  S),  and  a  similar 
solution  by  Frederick^^’  refer  specifically  to  underwater  shock.  It  is  proposed  now  to  use  the 
results  of  these  analyses,  particularly  that  by  Hudson,  to  predict  the  collapse  of  plates  due  to  blast 
effects  in  air.  These  predictions  will  then  be  subjected  to  experimental  checks  under  the  test  pro¬ 
gram  at  present  being  formulated.  If  it  is  found  that  the  discrepancies  between  theoretical  and 
experimental  results  are  sufficiently  serious,  a  fresh  theoretical  analysis  will  then  be  undertaken.  It 
is  believed  that  this  approach  will  lead  to  the  greatest  economy  of  effort  in  the  overall  program. 


Z.  Z.  Z.  1  Summary  of  the  Hudson  Approach 

The  problem  attacked  by  Hudson  is  the  damage  done  by  underwater  explosions  to  thin  metal 
circular  diaphragms,  air  backed,  and  held  rigijily  at  the  peripheries.  The  material  of  the  diaphragm 
is  supposed  to  be  rigid- perfectly  plastic.  Although  Hudson  has  also  given  an  approximate  solution  for 
a  work-hardening  material,  the  primary  difference  noted  is  that  in  the  work- hardening  solution  the 
apex  of  the  cone  (the  final  shape  assumed  by  the  diaphragm)  is  rounded  off,  whereas  in  the  perfectly 
plastic  material  the  apex  is  a  point.  For  the  present  purposes,  it  appears  sufficient  to  fix  attention 
on  the  solution  for  the  rigid- perfectly  plastic  material. 

The  shape  assumed  by  the  diaphragm  at  any  time  t  <  ts,  where  tg  is  the  total  time  for 
deflection  or  "swing  time,  "  is  that  of  a  truncated  cone  as  shown  in  Figure  Z.  3.  The  flat  central 
region  travels  with  a  constant  speed  v  normal  to  its  plane,  decreasing  in  radius  with  its  motion,  the 
diaphragm  finally  assuming  the  shape  of  a  cone  with  a  zero  thickness  at  the  apex. 

While  an  explicit  form  of  the  solution  of  the  equations  of  motion  is  not  given  even  for  the 
highly  idealized  model  proposed  by  Hudson,  certain  special  cases  can  be  solved  and  information 
gleaned  regarding  the  following: 

(1)  The  radius  R  of  the  bending  wave  (see  Fig.  Z.  3)  as  a  function  of  time 

(Z)  The  diaphragm  profile  at  each  instant 

(3)  The  thickness  distribution 

(4)  Dis  pi  a  cement -time  curves  of  particles  in  diaphragm 

(5)  Stress  and  strain  distributions 

(6)  The  center  deflection  as  a  function  of  v 

(7)  The  total  time  for  the  deformation  to  take  place 


The  solutions  given  by  Hudson  are  summarized  below. 
Z.Z.Z.Z  An  Elementary  Approximation  to  the  Solution 


When  the  acceleration  of  the  material  in  the  flat  central  portion  is  negligibly  small  through¬ 
out  the  motion,  i.e.,  U  is  a  constant,  the  solution  is  as  follows 


U  =  v^/Zc 

(a) 

R  =  a  -  ct 

(b) 

/R\ 

(c) 

Z  =  -  (a  -  R) 

(d) 

c 


(Z.83) 


The  thickness  distribution  in  the  deformed  diaphragm  given  by  (Z.  83  c)  shows  a  dimpling 
tendency  at  the  center.  In  fact,  at  the  last  moment,  the  thickness  becomes  zero  at  the  center.  How¬ 
ever,  as  the  variation  of  the  thickness  with  R  is  so  rapid  near  R  =  0,  the  tiny  pinhole  may  not  be  very 
apparent  experimentally. 
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FIGURE  2.  3.  A  SKETCH  OF  THE  DIAPHRAGM  CONFIGURATIONS 
AT  TWO  SUCCESSIVE  INSTANTS,  t  AND  t+dt.  THE  VARIOUS 
QUANTITIES  WHICH  ENTER  INTO  THE  SEVERAL 
GEOMETRIC  AND  KINEMATIC  RELATIONS 
ARE  PORTRAYED  HERE 
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Equation  (2.83  d)  shows  that  the  diaphragm  assumes  a  conical  shape  whose  center  deflec¬ 
tion  is  proportional  to  the  initial  velocity  v. 


The  total  time  of  deflection  tg  is  given  by 
tg  =  a/c 

2.  2.  2.  3  The  Exact  Solution  for  a  Rigid- Perfectly  Plastic  Material 

i(xi  -  l)/(xi-X2)  ■|(x2- 1)/(X2- xj) 

(a) 

\  (c2  .  xj  /  \  -  X2  / 

2  .  -2(xi)/(xi  -  X2) 


(2.84) 


f  =  P(hit) 


-  =  ^(i^)  = 
a 


(S) 


r* 


([x2  _  xj  )  (n2  .  X2) 


X 


P((i)dti 


m) 


-2(x2)/(x2-  xj) 


(b) 


(c) 


(2.85) 


where  xj,  X2  are  the  smaller  and  larger  (positive)  roots  respectively  of  the  quadratic  equation 
x^  -  (26  -  3)  X  +  2  =  0 

It  is  to  be  noted  that  the  initial  conditions  are 
P=q=l,  p=K  when  4=0 

Equation  (2.  85  c)  is  not  integrable  in  terms  of  elementary  functions. 

For  positive  values  of  6,  which  are  the  only  ones  of  physical  interest,  the  roots  of  (2.  86) 
occur  in  pairs  of  positive  values  for  ^  +  ’'Iz  ^  6  ^  oo  and  as  pairs  of  negative  values  for  0^6 
(a  very  small  range).  For  intermediate  values  of  6,  the  roots  are  complex.  Hudson  confine s  himself 
to  a  consideration  of  values  of  6  >  |  +  '77.  At  the  lower  limit  of  this  range,  corresponding  to  very 
large  v,  or  small  cr,  values,  we  have  xj  =  X2  =  *72  •  As  6  becomes  very  large,  n/2  >  xj  —  0  and 

>72  <  X2  —  CD. 

A  case  of  considerable  interest  arises  when  coincides  with  one  of  the  roots  of  (2.  86). 

From  (2.  85),  it  is  seen  that  the  only  possibility  is  for  itself  to  be  constant  and  equal  to 
.  It  can  then  be  shown  that  either  0<p=K<l.  orp=-K<-l.  The  first  case  only  is  admissible, 
since  the  second  precludes  p  =  k  initially#  Thus,  we  have 


0<p2=,2  =  ,j=  (6-|)  -  {(S-I) 


and 


6  >  3 
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The  solution  is  then  as  follows 


1  - 

R  =  a  -  V,  vt  =  a  - 

Zk 

1  - 

1-JL-  z 

Zk 

(a) 

1 

li 

X 

(b)  ^ 

U  =  KV 

(c) 

,  ,  K**  +  3(c2  +  2 

>  3 

(d) 

For  K  =  1,  6  =  3  and  v  =  (Z/’\[3)c,  and  the  deformed  diaphragm  is  cylindrical,  i.e.,  it  ruptures  com¬ 
pletely. 


For  other  allowed  values  of  k  (=|ji),  this  solution  is  quite  similar  to  the  elementary  approx¬ 
imation  presented  earlier,  even  quantitatively,  provided  6  is  large.  For  it  follows  then 

K  «  =  v/2c«l,  (6»  3) 


and  (2.  87)  reduces  to  (2.  83). 


Another  case  of  interest  arises  if  k  has  any  desired  positive  value  within  reason  and  6  is 
very  large.  Then  xj  is  small  and  x^  is  large.  Assume  that  K^/x2  «  1|  as  is  (jl^/x^  (since  varies 
from  fc^  to  xj).  Under  these  simplifying  conditions,  it  is  possible  to  derive  an  approximate  integral 
of  (2.  85  c).  It  is  then  found  that 


X]^  «  1/6,  X2  ^  26  »  1 


Then  providing  6  is  large  enough 


—  =  —  =  2Uv/(v2  -  u2)  [  1 

-  (R/a)] 

(a) 

a  a 

/R  \ 

(b) 

h  \  a  / 

(c) 

(2.  88) 


It  may  be  shown  that  this  approximation  is  uniform  over  the  range  of  U/v  from  K  to  v/2c,  and  the 
range  of  R/a  from  1  to  0;  but  the  nearer  k  is  to  1  the  larger  6  must  be.  If  K  =  0,  which  is  the  case  of 
a  material  which  has  strictly  no  elastic  strain  range,  the  deflected  shape  of  the  diaphragm  described 
by  (2.  88  a)  is  conical  near  the  center  with  a  central  deflection  the  same  as  that  given  by  the  elemen¬ 
tary  theory  under  2.  2.  2.  2. 


2.  2.  2.  4.  General  Remarks 


It  may  be  noted  that  the  modulus  of  elasticity  E  enters  into  the  definition  of  k  .  This  is 
because,  although  the  elastic  phase  of  the  material  is  neglected  in  the  specification  of  the  deforma¬ 
tions,  an  estimate  has  been  made  of  the  initial  conditions  at  the  end  of  the  elastic  range.  The  time 
at  which  elastic  phase  is  completed  is  of  the  order  of  magnitude  of  3%  of  tg  for  some  steels.  Thus, 
it  is  generally  immaterial  whether  time  is  reckoned  from  the  end  of  the  elastic  strain  range  or 
before  it. 


Whether  the  above  analysis,  developed  with  specific  reference  to  diaphragms,  would  apply 
to  plates  of  the  proportions  of  interest  in  the  present  research  program  can  only  be  determined  by 
careful  experimentation  and  photographic  records  of  the  deformation. 
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Some  additional  results  which  may  be  of  value  in  checking  the  theory  are 


C-J.  =  0-0  =  0- 

Cj.  =  eg  =  log  {r/r^) 
ejj  =  log  H/h 


1 


(2.89) 


31 


REFERENCES 


2.  1  Fung,  Y.  C.  and  Sechler,  E.  E.  ,  "Instability  of  Thin  Elastic  Shells,  "  Proc.  ,  First 

Symposium  of  Naval  Structural  Mechanics,  Pergamon  Press,  New  York,  19^0. 

2  Th.  von  Kerman,  "Encyklopadie  der  Mathematischen  Wissenschaften,  "  vol  IV^, 

p  349,  1910. 

2.3  G.  Hermann,  "Influence  of  Large  Amplitudes  on  the  Flexural  Motions  of  Elastic  Plates," 
NACA  TN  3578. 

2.4  Hu-Nan  Chu  and  G,  Hermann,  "Influence  of  Large  Amplitudes  on  Free  Flexural  Vibra¬ 
tions  of  Rectangular  Plates."  Journal  of  Applied  Mechanics,  p  532,  December  1956. 

2.  5  H.  M.  Berger,  "A  New  Approach  to  the  Analyses  of  Large  Deflections  of  Plates.  " 

Journal  of  Applied  Mechanics,  p  465,  December  1955. 

2.6  W.  A.  Nash  and  J.  R.  Modeer,  "Certain  Approximate  Analyses  of  the  Nonlinear 
Behavior  of  Plates  and  Shallow  Shells."  Proceedings  of  the  Symposium  on  the  Theory 
of  Thin  Elastic  Shells,  Delft,  August  1959,  Interscience  Publishers,  Inc.  ,  New  York, 
I960. 

2.7  N.  W.  McLachlan,  "Bessel  Functions  for  Engineers,  Edition,  Oxford,  1955. 

2.8  G.  E.  Hudson,  "A  Theory  of  the  Dynamic  Plastic  Deformation  of  a  Thin  Diaphragm,  " 
Journal  of  Applied  Physics,  Vol.  22,  No.  1,  Jan.  1951,  p  1. 

2.  9  Daniel  Frederick,  "A  Simplified  Analysis  of  Circular  Membrane  Subjected  to  an 

Impulsive  Loading  Producing  Large  Plastic  Deformations,  "  4th  Midwestern  Conference 
on  Solid  Mechanics,  1959,  p  18. 


32 


3.  NUMERICAL  ANALYSIS  -  CYLINDRICAL  SHELLS  (LINEAR- ELASTIC,  SMALL  DEFLECTIONS) 


The  computations  utilized  the  following  geometries  and  materials  for  the  cylindrical  models 
used  in  the  experimental  investigation 


a  =  6" 


V  =  1/3 


L  =  36" 
h  =  .036" 


E  =  2.  9  X  10  (psi) 

.-4 


p  =  7.  36  X  10" 


/ib-sec*^^ 

V  in4  / 


The  analyses  began  with  determining  the  free  vibration  circular  frequency  utilizing  the  values 
obtained  for  the  frequency  factor  (4j)  in  Equation  (I.  13)  for  the  appropriate  values  of  k  [  see  Eqs. 

(2.  9),  (2.  16)  or  (2.  17)1  in  the  relationship 

CO  =(12.34,  )^^^X  10^ 

mn 

The  variation  of  the  circular  frequency  with  the  end-closure  stiffness  (R)  for  various  values  of  m  and 
n  are  shown  in  Figure  3.  20. 

The  computations  in  all  cases  were  directed  toward  obtaining  pertinent  radial  displace¬ 
ments  from  the  relation 

w(x,  4),t)=  5^  <t>)  q^^(t) 

Values  for  q  (t)  in  Equation  (3.  1)  were  acquired  from  Equation  (I.  27)  I  or  in  finite  difference  form, 
Eq.  (1.28)]  ,  ^^lich  entailed  obtaining  the  generalized  force  and  generalized  mass. 

3.  1  Simply  Supported  End- Closures 

For  the  generalized  force  [  see  Eq.  (I.  27);  Appendix  l]  ,  taking  the  loading  as  being  uni¬ 
formly  distributed  along  the  length  of  the  shell,  we  have 


Qmn(t)  =  " 


f  [?(<!>,  t)  J  Wj^n^x,  <t>)dx1  d4> 

^  •'o  -* 


(3.  3) 


From  Equation  (2.  15) 


2tt 


2  EaO 


Qmn(t)  = 


mn'*''  nir 


/p(4>. 


t)  cos  m(t)d<t)  (n  =  1 ,  3,  5,  .  .  .  ) 


(3,4) 


Taking  into  consideration  the  loading  symmetry  and  using  cylinder  surface  increments  of  ^  22.  5  ,  ^ 

Equation  (3,4)  becomes 


4LaC. 


mn'  '  mr 


-  11.25 

/ 


P^Ct)  cos  m<t)d<j>  + 

/ 


1.1 


4>i 


+  22.  5 


►180 


168. 75 


Pl8o(^) 


] 


^4,^+  11.25  cosm4>d<|> 


(3.5) 


33 


The  values  of  P±(t)  for  time  increments  of  0.  1  msec  were  obtained  from  load  profiles  similar  to 
those  shown  in  Figure  I.  10,  Appendix  I.  The  time  variation  of  acquired  from  Equation 

(3.  5)  for  Pj  =  33  psi  and  t^  =  1  msec,  is  shown  in  Figure  3.  1. 

From  Equation  (Z.  15)  [  and  Eqs.  (Z.  11)]  ,  the  expression  for  the  generalized  mass  [  see 
Eq.  (I.  Z7),  Appendix  l]  becomes 


M 


mn 


LihpTraC^ 

Z 


(3.6) 


The  frequencies  and  the  corresponding  values  of  the  terms  within  the  brackets  of  Equation  (3.  6)  are 
given  in  Table  3.  1. 


In  each  of  the  computations,  only  the  lowest  of  three,  real  frequencies  obtained  from 
Equation  (I.  13)  was  used.  This  frequency  was  associated  with  the  response  where  the  radial  dis¬ 
placement  (w)  was  predominant.  An  indication  of  the  contributions  made  by  the  longitudinal  (u)  and 
tangential  (v)  displacements  may  be  found  in  Table  3.  1  where  sample  values  of  K^/K^  (reflecting  the 
tangential  displacements)  and  K^mr/K^L  (reflecting  the  longitudinal  displacements)  are  given.  Here, 
it  is  seen  that,  although  these  values  are  greater  than  unity,  they  are  at  the  same  time  associated 
with  much  higher  frequencies.  Since  the  value  of  the  generalized  mass  is  increased  [  Eq.  (3.  6)]  and 
the  larger  frequencies  dictate  smaller  values  of  At  in  Equation  (I.  Z8),  the  net  result  is  a  drastic 
decrease  in  the  values  of  This  clearly  shows  that  the  contribution  of  the  longitudinal  and 

tangential  displacements  is  quite  negligible  when  compared  with  that  of  the  radial  displacement. 


The  ratio  of  the  generalized  force  to  the  generalized  mass  appearing  in  Equations  (I.  Z7) 
and  (I.  Z8)  is,  in  accordance  with  Equations  (3.  5)  and  (3.  6) 


Q  (t) 

mn^  ' 

M 


nhpTT 


P^{t)  cos  m<|)d(|) 


(m  =  2,  3, .  .  .  ) 
(n  =  1,3,5,...) 


(3.7) 


where,  for  the  particular  shell  properties  previously  listed 


4 


nhpTT 


Z 


1.  53  X  10*^’ 
n 


From  Equation  (I.  Z8),  the  variations  of  qj^j^(t)  with  time  are  obtained  for  particular  values  of  Pj  and 
tf.  These  are  shown  in  Figures  3.  Z  through  3.  13.  As  a  final  step,  Equation  (3.  Z)  is  utilized  to 
obtain  the  value  of  the  radial  displacement  for  each  point  on  the  cylinder’s  surface  as  defined  by  x 
and  <|>.  In  Figures  3.  Z  through  3.  13  are  shown  the  variation  of  w(x,  (\>,t)  with  time  for  the  point  at 
x  =  L/Z  and  <|>  =  0°.  (See  adso  Table  3.  Z) 

3.  Z  Fixed  End-Closures 


For  cylinders  with  fixed  boundary  conditions,  the  expression  for  the  generalized  force 
[see  Eqs.  (Z.  13)  and  (3.  3)]  for  a  load  uniformly  distributed  along  the  length  of  the  shell  is 


i(t) 


LaC 


P(k) 


P{(j),  t)  cos  m<|)d<|) 


(3.8) 


where 

P(k)  =  r^(cosh  K  +  cos  fc  -  Z)  +  (sin  k  -  sinh  /c) 
F’l  =  sinh  /c  -  cosh  #c  +  sin  k  +  cos  #c 


34 


c/  j 
JtlZ. 
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VARIATION  OF  LOADING  FUNCTION  WITH  TIME 


TABLE  3.  1.  FREQUENCIES  AND  RELATED  QUANTITIES  FOR  CYLINDRICAL 
SHELLS  WITH  SIMPLY  SUPPORTED  END-CLOSURES 


CO 

< 

M  8 


vO 


00  nO 


<N3 


t: 

d 

cn 


rq  rq 
m  — «  o 

rs3 

vO 


<N3 

(M 


O 

00  ^ 
nO  00  *-1 


o 


lO  *—«  lO  •— I 

00 


c 

CO 


(Ml  —I 


CO  CO 
—t  r- 


00 

r- 


o  o  m 
o  LO  r- 


00  vO 

-"t  o 

(vi 

o 

(M 


CO  CO 
O  CO  CO 


ir> 


in 

(VJ 


o 

O 

o 

O 

o 

o 

o 

o 

o 

00 

(M 

CO 

in 

00 

o 

CO 

00 

o 

CO 

vO 

CO 

00 

(M 

vD 

r- 

o 

CO 

u 

o 

o 

O 

o 

O 

o 

o 

O 

o 

0) 

CO 

(M 

CO 

o 

•-H 

CO 

in 

o 

in 

(D 

O 

r- 

o 

•V 

.V 

00 

ro 

o 

o 

•-H 

o 

CO 

00 

in 

(d 

r- 

■-H 

sO 

(M 

u 

r-H 

61 


X 

in 

O 

u 

0) 


o 

o 

o 


l-H 

X 

X 

■H 

•H 

73 

73 

d 

d 

a> 

(0 

a 

a, 

a 

a. 

< 

< 

fM 

o 

00 

•-H 

(M 

»-4 

CO 

6* 

a" 

6* 

W 

W 

w 

(U 

0) 

<0 

d) 

(U 

(D 

w 

w 

t— 1 

(M 

CO 

36 


S/MPL.'V  ^  SL/PPO/2  7~£.£^  £A>/G  <ZLGSUP.£S 

*  /7  pS/ 

^  s  0.7&  /77^<SC. 

_  £xC^l.UJ9tA^G  /A/ 

_ /f^cLiJGfh/3  u4>^  Mmr> 


-tv  (i)/^.i/.  ,  ^ 


HI 


.3  ,4  .5  .6  "7  /•<=»  >-f  ^*2 

TlK4e^  (/JO-SCC^) 

FIGURE  3.  2.  VARIATION  OF  AND  RADIAL  DISPLACEMENT  w(t) 

AT  X  =  L/2,  <t>  =  0“  FOR  CYLINDRICAL  SHELL 
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FIGURE  3,  3.  VARIATION  OF  qjnn(t)  AND  RADIAL  DISPLACEMENT  w(t) 
AT  X  =  L/2.  ({)  -  0“  FOR  CYLINDRICAL  SHELL 


FIGURE  3.4.  VARIATION  OF  qinn(t)  AND  RADIAL  DISPLACEMENT  w(t) 
AT  X  L/2,  <})  =  0“  FOR  CYLINDRICAL  SHELL 


FIGURE  3.  5.  VARIATION  OF  RADIAL  DISPLACEMENT  w(t) 

AT  X  =  L/2,  <j)  =  0“  FOR  CYLINDRICAL  SHELL 
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FIGURE  3.  6.  VARIATION  OF  q^n^t)  AND  RADIAL  DISPLACEMENT  w{t) 
AT  X  L/3,  <|>  ^  O'  FOR  CYLINDRICAL  SHELL 
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FIGURE  3.7  VARIATION  OF  qmn(t)  ^.ND  RADIAL  DISPLACEMENT  w(t) 
AT  X  L/Z,  4)  ^  0“  FOR  CYLINDRICAL  SHELL 
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FIGURE  3.  10:  VARIATION  OF  qmn(t)  AND  RADIAL  DISPLACEMENT 
AT  X  =  L/2,  <})  =  0"  FOR  CYLINDRICAL  SHELL 
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FIGURE  3.  II.  VARIATION  OF  qinii(t)  AND  RADIAL  DISPLACEMENT  w(t) 
AT  X  =  L/2,  <|>  =  0“  FOR  CYLINDRICAL  SHELL 
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FIGURE  3.  12.  VARIATION  OF  qmn(t)  AND  RADIAL  DISPLACEMENT  w(t) 
AT  X  =  L/2,  <|>  -  0"  FOR  CYLINDRICAL  SHELL 
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FIGURE  3.  13.  VARIATION  OF  qmn(t)  AND  RADIAL  DISPLACEMENT  w(t) 
AT  X  =  L/2,  6  =  0'’  FOR  CYLINDRICAL  SHELL 


TABLE  3.  2.  COMPUTED  MAXIMUM  RADIAL  DISPLACEMENTS  (COLUMNS  A)  AND 
CORRESPONDING  TIME  OF  OCCURENCE  (COLUMNS  B)  AT  x  =  L/2,  c|>  =  0° 
FOR  CYLINDERS  WITH  SIMPLY-SUPPORTED  END- CLOSURES 


Columns  A  _ Cdlumns  B 


Pi 

t 

+ 

(1) 

(2) 

(1) 

(2) 

I 

.2 

(3) 

(psi) 

sec) 

(inches) 

(inches) 

(m! 

sec) 

jmj 

sec) 

lb  -in^ 

-msec 

17 

0. 

75 

.  042 

.  035 

0. 

,  9 

0. 

85 

6. 

3 

17 

2. 

5 

.  150 

.  085 

L 

.  4 

1. 

4 

21. 

2 

23 

1. 

0 

.  078 

.  065 

0, 

.  9 

0. 

9 

11. 

5 

25 

2. 

0 

.  114 

.  092 

1, 

.  3 

1. 

3 

25 

25 

3. 

0 

.  152 

.  124 

1 

.  4 

1. 

4 

37. 

5 

33 

1. 

0 

.  084 

.067 

1 

.  0 

1. 

0 

16. 

5 

35 

2. 

0 

.  164 

.  128 

1 

.  4 

1. 

3 

35. 

0 

37 

L 

5 

.  132 

.  107 

1 

.  3 

L 

3 

27. 

7 

37 

3. 

0 

.  212 

.  174 

1 

.  3 

1. 

3 

55. 

5 

60 

1. 

0 

.  167 

.  133 

1 

.  0 

1. 

0 

30 

60 

2. 

0 

.  284 

.  184 

1 

.  4 

1. 

,  4 

60 

60 

3. 

,  0 

.  325 

.26 

1 

.  3 

1. 

2 

90 

(1)  Neglecting  contributions  of  u  &  v  in  Mj^ji 

(2)  Including  contributions  of  u  &:  v  in  Mj-Qn 

(3)  Pi  X  4/2 
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^2 

=  sin  K  +  sinh  k 

^3 

=  cosh  K  -  cos  K 

K 

=  n*Tr 

n' 

=  l/2(2n+  1) 

(n  =  1,2,3,.. 

.) 

It  is  of  interest  to  note  that  the  integral  term  in  Equation  (3.  8)  is  identical  to  that  appearing  in  Equa¬ 
tion  (3.4).  Thus,  the  values  of  the  generalized  force  for  the  simple  and  fixed  support  conditions  are 
proportional.  The  values  of  the  terms  in  Equation  (3.  8)  are  given  in  Table  3.  3. 

For  the  generalized  mass,  in  a  first  order  approximation,  the  longitudinal  and  tangential 
components  of  the  displacement  were  neglected  [  see  Eq.  (I.  Z7),  Appendix  l]  ,  so  that 


M 


=  pah  J  J  f 


or 


M 


mn 


LhpairC^ 

(r”) 


(3.9) 


whe  re 


a^(K)  =  sinh  k  cosh  k  +  k  ~  cosh  k  sin 

-  ^2^3  cosh  2k  -  2  sinh  k  sin  k  -  --  cos  2k ^ 

2/1  \ 

+  r 2  I  sinh  K  cosh  k  -  cosh  k  sin  k  j 

The  appropriate  values  of  a^(K:)  are  also  given  in  Table  3.  3 

From  Equations  (3.  8)  and  (3.  9),  we  acquire  the  ratio  necessary  for  the  solution  of  Equa¬ 
tion  (I.  28) 

^  P(<c)  V  r  (m  =  2,3,...) 

-  phirC  ai(K)  J  P4,(t)  cos  m4.d<|>  (n  =  1 , 3,  5.  .  .  .  )  (3.10) 

<t> 


where  the  summation  is  identical  to  the  terms  within  the  brackets  in  Equation  (3.  5).  The  variation 
'lmn(l)  time,  obtained  by  utilizing  Equation  (3.  10)  in  Equation  (I.  28)  are  shown  in  Figures 

3.  14  through  3.  19.  The  variation  of  the  radial  displacement  with  time  (also  shown  in  Figs.  3.  14 
3.  19)  are  obtained  from  Equation  (3.  2),  or,  for  the  fixed  boundary  conditions 


w(x,<|).  t)  =  ^  ^  Fj  (cos^  -  cosh^)  +  (  sinh  ^  -  sin 

Table  3.4  contains  the  maximum  radial  displacement  and  corresponding  times. 


cos  m((>q(t) 


(3.12) 


In  order  to  determine  the  contributions  of  the  longitudinal  (u)  and  tangential  (v)  components 
of  displacements,  we  use  the  expression  for  generalized  mass  appearing  in  Equation  (I.  27)  and  from 
Equations  (2.  11)  and  (2.  12)  obtain 
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TABLE  3.  3.  FREQUENCIES  AND  RELATED  QUANTITIES  FOR 
CYLINDRICAL  SHELLS  WITH  FIXED  END-CLOSURES 


^mn 

mn 

m 

n 

n' 

(rad/ sec) 

cycle/ sec 

1 

1 

1.  5 

50,000 

7,960 

220 

1 

2 

2.5(3) 

1.  0 

1 

3 

3.  5 

56, 310 

8,960 

120,000 

2 

1 

1.  5 

4,  640 

740 

220 

2 

2 

2.5(3) 

13,  660 

2,  180 

1.  0 

2 

3 

3.  5 

79, 730 

12, 700 

120,000 

3 

1 

1.  5 

2,  040 

325 

220 

3 

2 

2.5(3) 

6,  010 

960 

1.  0 

3 

3 

3.  5 

12,050 

1,920 

120, 000 

25, 600 

11,170 

9684 

.4525 

25  X  10^ 

0 

,  8266 

.  6167 

9.  7  X 10^^ 

264.  8  X  10^^ 

,  6338 

.  6172 

25, 600 

25  X  10° 

9.  7  X 10^^ 

11,170 

0 

264.  8  X  10^^ 

,  5022 

.  2007 

25, 600 

25  X  10^ 

9.  7  X 10^^ 

11,170 

0 

264.8  X 10^^ 

.  3335 

.  0874 

(1)  See  Eq.  3.  8 

(2)  See  Eq.  3.  9 

(3)  Anti  symmetrical 

(4)  See  Eq.  3.  13 

(5)  aiM 


(5) 

51895 

32506 

28530 
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FIGURE  3.  15.  VARIATION  OF  AND  RADIAL  DISPLACEMENT  w(t) 

AT  X  =  L/2.  <}>  =  0“  FOR  CYLINDRICAL  SHELL 
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FIGURE  3.  18.  VARIATION  OF  qjnn(t)  AND  RADIAL  DISPLACEMENT  w{t) 
AT  X  =  L/2,  A  =  0^  FOR  CYLINDRICAL  SHELL 
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FIGURE  3.  19.  VARIATION  OF  qj^j^(t)  AND  RADIAL  DISPLACEMENT  w(t) 

AT  X  =  L/2,  (})  =  0^  FOR  CYLINDRICAL  SHELL 
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TABLE  3.  4.  COMPUTED  MAXIMUM  RADIAL  DISPLACEMENTS  (COLUMNS  A)  AND 
CORRESPONDING  TIME  OF  OCCURENCE  (COLUMNS  B)  AT  x  =  L/2,  (j)  =  0“ 
FOR  CYLINDERS  WITH  FIXED  END-CLOSURES 


Columns  A  _ Columns  B 


^I 

(1) 

(2) 

(1) 

(2) 

I  (3) 

(£Si) 

(msec) 

(inches) 

(inches) 

(msec) 

(msec) 

(lb  -in  -msec) 

17 

0.  75 

.  030 

.  025 

0.49 

0.  47 

6.  3 

17 

2.  5 

.  039 

.  035 

0.  54 

0.  54 

21.  2 

33 

1.  0 

.  051 

.  040 

0.  46 

0.  46 

16.  5 

35 

2.  0 

.  068 

.  057 

0.  60 

0.  53 

17.  5 

37 

1.  5 

.  049 

.  040 

0.  43 

0.  48 

27.  7 

60 

2.  0 

.  093 

.  076 

0.  53 

0.  53 

60 

(1)  Neglecting  contributions  of  u  and  v  in  M^n 

(2)  Including  contributions  of  u  and  v  in  M^^^ 

(3)  (Pi)(t+) 

2 
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where 


a2^(ic)  =  ^  3  (  sinh  k  cosh  k  +  cosh  k  sin  k)-  j  (cosh  Zk  +  cos  Zk) 

2/1  \ 

+  r  ^  I  —  sinh  K  cosh  fc  +  /c  -  cosh  k  sin  k  j 


(3.  13) 


+  r  ^  ^  —  sinh  K  cosh  fc  +  /c  -  cosh  k  sin  k  J 
Accordingly,  Equation  (3.  10)  is  replaced  by 

Q  p„  P(k)  V  /  P^(t)  cos  m((>d(|> 

^mn^  '  -  ^  4) _ 

^mn  "  r  /T. 


p 

/  r,  \  ^“1 

/  x  2 

f  \ 

fK  K  \ 

aj(K)  1  + 

{ _ 

+  viqx  y 

(3.  14) 


The  corresponding  values  for  q^n^^)  obtained  from  Equation  (I.  28)  and  w(t)  at  x  =  E/2  and  (p  =  0°  are 
given  in  Figures  3.  14  through  3.  19. 

3.  3  Elastic  End-Closures 

For  cylinders  with  flexible,  flat-plate,  end-closures,  we  take  as  a  first  order  approxima¬ 
tion  the  following  relation  for  the  moment 


M  d<t)  ^  ^  <|) 


(3.  15) 


D  dr  r 

where  the  bar  denotes  properties  associated  with  the  plate.  The  solution  of  Equation  (3.  15)  for  <\>  is 

^  _  M  T 
9  -  —  2 

D  v+  1 

Setting  the  value  of  (|)  at  r  =  a  equal  to  one  radian,  the  value  of  stiffness  Sis  found  by 
o  5(1  +15) 

S  =  -  (for  <|)  I  =1  radian) 

a  la 

If  the  flat  plate  has  the  same  thickness  as  the  cylinder  wall  and  is  made  of  the  same  material  as  the 
shell,  then  using  the  specific  dimensional  and  materials  values  previously  listed,  we  have  for  this 
particular  case 


and  thus 


S  =  28.  2 


F or  a  load  uniformly  distributed  along  the  length  of  the  shell,  the  expression  for  the  gener¬ 
alized  mass  is 
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LaC 

Qmn^t)  =  TTT-  P(k,R)  J  P(<t>.  t)  cos  m<{>d4) 


(3.  16) 


whe  re 


P(»c,R)  =  2  [  sin  K(cosh  k  -  1)  +  sinh  fc(l  -  cos  k)  +  R/cCcosh  k  cos  k  -  1) 

+  R^[  sin  ic(cosh  fc  -  1)  -  sinh  #c(l  -  cos  k)]  J 
=  2k^  sin  K  +  R(k  +  R)(sinh  k  -  cosh  k  +  sin  k  +  cos  k) 


=  Zk^  sin  fc  +  R(fc  +  R)  r^’  [  See  Eq.  (3.  8)] 

r  2  =  sin  K  +  sinh  k 
r  ^  =  cosh  K  -  cos  K 

and  K  are  the  values  commensurate  with  the  value  of  R  obtained  from  Figure  2.  1.  These  values  of  k 
and  the  circular  frequencies  obtained  from  Figure  3.  20  are  given  in  Table  3.  5. 

For  the  generalized  mass,  neglecting  the  contributions  of  the  longitudinal  and  transverse 
components  of  displacements,  we  have 


^  Lh£^  ajU.R) 


(3.17) 


where 


aj(/c,  R)  =  2K:^r2(sinh  k  -  sin  k)  +  Zk^I  3  sinh  fc  sin  /c(cosh  k  sin  k  ~  sinh  k  cos  k) 

3  2 

-  RTj  (sinh  k  -  sin  k)]  +  k'^R[  4  sinh  k  sin  K(cosh  k  cos  k  -  1)  -  RT 2I 

+  k"r"[  (cosh  K  cos  K  -  l){sinh  k  cos  k  •*  cosh  k  sin  /c)  +  8  sinh  k  sin  k  (sinh  fc  cos  k 

+  cosh  K  sin  k)  +  Zr^P^  ]  +  /cR^l  Rr^+  4  (cosh  k  cos  k  -  1)^1 

+  3Ri[cosh  K  cos  K  -  l)(cosh  fc  sin  k  +  sinh  k  cos  k) 

From  Equations  (3.  16)  and  (3.  17),  the  necessary  ratio  for  the  solution  of  Equation  (1.  28)  is  as 
follows 


phife  tI^)  Z  / 


(3.18) 


where  the  summation  is  identical  to  the  terms  within  the  brackets  in  Equation  (3.  5).  The  variation 
of  w(x,<(i,  t)  at  X  =  L/2,  and  <()  =  0®  with  time  is  shown  in  Figures  3.  21  through  3.  23. 
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TABLE  3.  5.  FREQUENCIES  AND  RELATED  QUANTITIES  FOR 
CYLINDERS  WITH  ELASTIC  END- CLOSURES 


m  n  K 

2  1  1.4117 

3  1  1.41tt 


C*} 

mn 

(rad/ sec) 
4,  000 
2,  100 


p(fc, 

12, 200 
12, 200 


56.  49  X  10^ 
56.  49  X  10^ 


(1)  See  Equation  (3.  16)  (for  R  =  8) 

(2)  See  Equation  (3.  17)(for  R  =  8) 
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FIGURE  3  22.  VARIATION  OF  RADIAL  DISPLACEMENT  *.0(1) 

AT  X  =;  L/2,  <|>  =  0®  FOR  CYLINDRICAL  SHELL 


FIGURE  3.  23.  VARIATION  OF  AND  RADIAL  DISPLACEMENT  «-(t) 

AT  X  =L/2,  (j)  ^0®  FOR  CYLINDRICAL  SHELL 

5^i 


4.  DISCUSSION 


The  utility  of  an  analytical  procedure  is  best  measured  in  terms  of  its  ability  to  yield 
meaningful  and  correlative  quantities  for  a  wide  range  of  variables  with  a  minimum  of  computational 
effort.  The  adequacy  of  these  computed  quantities  can  be  ascertained  only  in  a  direct  comparison 
with  reliable,  equivalent  experimental  evidence.  If  the  results  obtained  by  a  comparatively  simple 
analytical  method  agree  favorably  with  the  available  experimental  data  and  if  the  results  cover  the 
range  of  interest  for  the  problem  at  hand,  there  then  appears  to  be  little  reason  (at  least  from  a 
practical  viewpoint)  for  employing  another  more  rigorous  and  therefore  more  computationally  com¬ 
plex  solution. 

The  linear- elastic  method  for  cylindrical  shells  developed  in  Section  2.1.1  and  used  in  the 
numerical  analysis  in  Section  3  appears  to  yield  significant  results  (over  a  wide  range  of  peak  over¬ 
pressures  and  positive  phase  durations)  that  agree  favorably  with  the  experimental  data.  There  is, 
at  present,  no  apparent  or  justifiable  reason  for  suggesting  the  use  of  other  available  but  neverthe¬ 
less  more  complex  analytical  procedures  (such  as  that  developed  in  Section  2.1.2). 

One  method  of  compactly  presenting  the  numerical  results  obtained  in  the  preceding  sec¬ 
tion  is  to  chart  the  variation  of  the  maximum  radial  displacement  (in  this  case,  at  x  =  L/2,  and 
(j)  =  0®)  with  the  two  characteristics  which  distinguish  each  particular  blast  loading,  namely,  the 
peak  incident  overpressure  (Pj)  and  the  companion  positive  phase  duration  (t^).  In  Figure  4.  1,  we 
have  taken  the  impulse  (defined  as  Pj  X  t^  -r  2)  for  the  ordinate  axis  and  the  maximum  radial  deflec¬ 
tion  for  the  abscissa.  On  one  hand,  the  bands  for  each  type  of  cylinder  end-closure  arise  from 
taking  the  generalized  mass  as  a  function  of  only  the  radial  component  of  displacement  and  on  the 
other  including  all  three  displacement  components  (u,  v  and  w)  in  computing  Table  4.  1  are 

comparisons  of  high-speed  camera  data  obtained  for  specific  models  with  the  computational  results 
as  represented  in  Tables  3.  2  and  3.  4  or  Figure  4.  1. 

The  close  agreement  between  the  experimental  and  analytical  results  is  indeed  encour- 
aging  particularly  in  view  of  some  of  the  expected  deviations  in  the  peak  incident  overpressure  and 
positive  phase  duration  that  arise  from  minor  differences  in  charge  size  and  shape,  differences  in 
atmospheric  conditions  and  local  changes  at  ground  zero  resulting  from  repeated  shots.  Moreover, 
there  are  the  assumptions  used  in  developing  the  load  spectrums  and  the  purely  analytical  approxi¬ 
mations  (see  Appendix  I)  and  some  of  the  normal  uncertainties  associated  with  the  high  speed  photog* 
raphy  data  (see  Appendix  III)  and  the  models'  dimensional  and  materials  properties.  All  these  fac¬ 
tors  undoubtedly  contribute  to  the  development  of  a  range  of  displacement  values  rather  than  a 
precise  and  specific  value  for  each  given  set  of  conditions.  The  encouraging  aspect  is  that  the 
numerical  procedure  used  in  the  preceding  section,  which  is  comparatively  simple  and  straight¬ 
forward,  does  give  values  that  are  compatible  with  the  experimental  evidence  and,  therefore,  of 
direct  practical  importance. 
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HALF  OF  PRODUCT  OF  PEAK  INCIDENT  OVERPRESSURE  AND 


POSITIVE  PHASE  DURATION 
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TABLE  4.  1.  COMPARISON  OF  EXPERIMENTAL  (HIGH-SPEED  PHOTOGRAPHY) 
AND  COMPUTED  (LINEAR- ELASTIC,  SMALL  DEFLECTION  THEORY) 
RADIAL  DISPLACEMENTS  FOR  CYLINDRICAL' SHELLS 


Computed _  Experimental^*^^ 


Model 

No. 

Shot 

No. 

Boundary 

Restraint 

Pi 

(psi) 

(msec) 

(PjXt+)/2 
(lb  -  in^  -  msec) 

(5) 

(inches) 

(6) 

(msec) 

(5) 

(inches) 

(6) 

(msec) 

CA3 

D1 

Fixed 

15  (2) 

6.  3  (4) 

47 

.064 

-  .078  (7) 

.04 

0.7 

CA3 

D2 

II 

19  (2) 

5.  5  (4) 

41 

.060 

-  .072  (7) 

- 

.07 

1.5 

CA3 

D3 

n 

42  (2) 

3.8  (4) 

80 

.098 

-  .  118  (7) 

- 

.  12 

0.  8 

CB2 

G1 

s,  s. 

20  (1) 

2.4(3) 

24 

.084 

-  .116  (7) 

1.4  (10) 

.  12 

1.3 

CB2 

G2 

II 

16  (1) 

2.3  (3) 

19 

.072 

-  .095  (7) 

1.4  (10) 

.  11 

1.  5 

CB2 

G3 

II 

17  (1) 

2.5  (3) 

21. 5 

.085 

-  .  105  (8) 

1.4(8) 

.  11 

1.4 

CB2 

G4 

II 

17  (2) 

2.9  (3) 

24.  5 

.088 

-  .  118  (7) 

1.4  (10) 

.  13 

1.3 

CB2 

HI 

II 

25  (2) 

5.1  (3) 

64 

.  194 

-  .248  (7) 

- 

.  20 

2.0 

CB3 

H3 

II 

23  (1) 

4.  5  (3) 

52 

.  164 

.  .  216  (7) 

- 

.  15 

1.5 

CB3 

H4 

II 

53  (1) 

3.9  (3) 

100 

.284 

-  .  366  (8) 

1.4(8) 

.30 

2.  0 

(1)  From  Table  III.  1,  Appendix  III 

(2)  From  Figure  III,  3,  Appendix  III 

(3)  From  Table  III.  2,  Appendix  III 

(4)  Computed 

(5)  Maximum  (inward)  radial  displacement  atx=L/2,  (t)=0® 

(6)  Time  of  occurrence  of  displacement 

(7)  From  Figure  4.  1 

(8)  From  Table  3.  2  or  computed 

(9)  From  Figures  III.  4  through  III.  13,  Appendix  III 

(10)  Estimates 
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u,  V,  w 


X,  <t),  Z 


METHODS  OF  ANALYSIS  (LINEAR- ELASTIC, 

SMALL  DEFLECTION  THEORY) 

LIST  OF  SYMBOLS 

Constant  coefficients  (maximum  amplitudes  of  component  vibrations) 
Eh^ 

12(1  -  v^) 

Elastic  modulus 
pa^(l  -  v^) 


E 

Eh 


Cylinder  length 
Bending  moment  of  section 
Generalized  mass 
Stress  resultant 
Generalized  force 
Kinetic  energy 
Potential  (strain)  energy 
Mean  cylinder  radius 
Cylinder  wall  thickness 

Subscript 


Number  of  circumferential  waves 

Generalized  coordinate 

acj) 

Time 

Longitudinal,  tangential  and  radial  components  of  displacements 
Axial,  tangential,  and  radial  coordinates 
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LIST  OF  SYMBOLS  (Cont*d) 


\/L 

Shear  strain 
Axial  strain 
Poisson's  ratio 


Circular  frequency 


I.  1  FREQUENCY  AND  DISPLACEMENT  EQUATIONS 


The  development  of  the  necessary  free-vibration  equations  follows  the  format  estab- 
lished  by  Yu.  (I-  The  solution  begins  with  the  following  equations  of  motion(  I-  2)  based  on  the 
assumption  that  -  0.  (See  Figs.  I.  1  and  I.  2  for  notation.  ) 


aN  aN. 

XX  ,  9X 

a  — ^ -  pahu  =  0 


dx 


a(j) 


9N  .  aN.  . 

X9  99 

a  +  pahv  =  0 


dx 


a  4> 


a  o  "i  r>  o  X  *1 


a  1  a  M^^ 


(i.i) 


Based  on  the  condition  that  a»z  (so  that  a  +  z  S  a)  and  the  following  equations  for  strains  as 
functions  of  displacements 


du 


a^w 


^x  = 


ax  -  ^  3x2 


1  dv  w  z  a^w 

"  a  a<t>  a  ■  ^2  g^2 


^X<t> 


dv  1  dv  ^  z  a  w 


dx  a  a<|)  a  dxdcl) 
we  obtain  the  following  differential  equations  for  the  displacements 


(1.2) 


4  1  a  w 

V^u  +  -  - = 


dxdy 


V  a^w  F  a^  fa  -  V  2v  aw 

ax^  'a^  at^  a(l-  V)  ax 


2F  a  u 
a^(l  -  v)  at^ 


].. 


(1.3) 


_4  2  +  V  a^w  1  a^w 

V  V  -  ~  2  “  ”  Q 

^  ax'^ay  ^  ay^ 


JL 

a^  atH'- 


■  -  V^v 

V 


dw 


2F  a  V 


a(l  -  v)  ay  a^(l  -  v)  at^ 


] 


=  0 


(1.4) 
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D  „8  1  - 

a^w  F 

8x4  - ^2 

^  3  -  V  2 

^  1  -  V 

2F  a^\  /  F 

’  a^(l  -  v)  8t^'  '  a} 

w 

°  V74  A 

4 

-  V  W  - 

2 

^ + ^i\ 

^  V  wj 

a^(l  -V) 

(I.  5) 


where 


y  =  acj) 


F  = 


pa^(i  -  vn 


D  = 


K  = 


12(1  -  vn 

Eh 


1  - 


The  space  and  time  dependent  forms  of  the  displacement  equations  are  taken  to  be 


ZXCL  * 

Aje  J  cos  m4)  sin  o^nt 

j 

ZXCl  * 

BjC  J  sin  mt|>  sin  «>mnt 

j 

’"mn<^’ ‘I’’ =  X  “mn* 

j 


y  (j  =  1.2,....8) 


(1.6) 

(1.7) 

(1.8) 


where 


J  L 


m  =  number  of  circumferential  waves 

(2m  =  number  of  circumferential  nodes) 


Substituting  Equations  (1.6),  (1.7),  and  (I.  8)  into  Equations  (1.3),  (1.  4),  and  (I.  5),  we  have 


^  [.  -  (^)  ]  [■-(¥)]} 
,  „  Z  [l  -  (2  *  V,  (i)  ]}  =  0 


(1.9) 
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[>  -  ( ^)  ]} 

aS  ajCje-jIm^  [>  ♦  -  (^)  ]  '  '  » 


(I.  10) 


Z  Cje^'^j  -  4'^n|^2  +  {3  -  v)m^  [l  -  (-J^r)]^' 

+  [l  -  (i?)  ]  +  [l  -  (■^)  ] 

+  (3  -v)km6  [l  -  (^)]  j-  (^  -  [^  -  {-^)] 

-  (1  -  v)(l  -  v^)(aaj)^j  =  0 


+  (3  -  v)km' 


-  (1  -  v)(l  -  v‘')(aa. 


(I.  11) 


where 


.  _  1  -  V  2  2 

^mn  E  P^mn^ 


Equations  (1-9),  (I.  10),  and  (I.  11)  are  simplified  by  assuming  that  the  shell  wall 
thickness  (a)  and  length  (L)  are  such  that 

C^)-'  « 

This  simplification  permits  Equations  (1.6),  (1.7),  and  (1.  8)  to  be  rewritten  as  follows 


Umn(x-  ‘t>.  t)  =  cos  mcf)  sin  cOj. 


=  tT"  ‘%i 


nt  ^ 


^tZcje^'^j 

j 


(j  =  l,2.3,4) 


(I.  12) 


Wmn(x,  4>.  t)  =  cos  me})  sin  ^ 
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where 


Kl  =  2‘1'mn  '  <3  -  v)  +  (3  -  v) 

=  m[m^{l  -  v)  - 
K3  =  +  (1  -  v)  m^] 

From  Equation  (I.  11),  we  also  obtain  the  following  frequency  equation 
2‘i'mn  -  ‘l'mn[3  +  (3  -  v)m^  +  2km^]  +  -  v)m^  +  (1  -  v) 

+  (3  -  v)km^]  -  (1  -  v)km®  -  (1  -  v)(l  -  v^)(aaj)^  =  0  (j  =  l,2,3,4)  (I- 13) 


The  solution  of  Equation  (I.  13)  for  (and  therefore  the  frequency  cOj^n)  depends 

on  obtaining  the  appropriate  values  of  aj(j  =  1,2,3,  4),  which  reflect  the  shell’s  end-closure 
boundary  condition.  For  the  fourth  degree  equation,  the  expected  values  of  Uj  would  be  of  the 
form 


K  K.  i-K 

'll  =  L’  "2  =  -  L’  “3  =  T  “4 


ifc 

TT 


where  i  =  sT-l  and  fc  is  a  real  number  containing  a  term  reflecting  the  axial  mode  shape. 

I.  2  Lagrangian  Equations  of  Motion 

The  development  of  the  Lagrangian  equations  of  n^otion  necessary  in  the  numerical 
analyses  for  the  response  of  cylindrical  shells  to  impulsive  loading  begins  with  the  general  dis- 


placement  equations 

u(x,  <t),  t)  = 

v{x,  <j),  t)  = 

Z  Vmn(’'’  <!>>*) 

► 

(I.  14) 

w(x,  <(),  t)  = 

^  ‘t*’  t) 

Utilizing  Equations  (I.  12),  these  may  be  rewritten  in  terms  of  generalized  coordinates 

u(x,  (|>,  t)  = 

Z  U^^{x,  4>)qnin(t) 

v(x,  4),  t)  = 

► 

(1.15) 

w(x,  4),  t)  = 

Z  Wjnn(x,  4^)qtnn(^)  j 
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The  expression  for  the  kinetic  energy  is 


L  ^Ztt  -  h/2 


T  =  i 


/  /  / 

0  0  -h/2 


[u(x,  <|),  t)]^  +  1  v{x,  (t>,  t)]  ^ +[w(x,  ((),  t-)]  ^ 


adxd(|)dz  (I.  16) 


Since  u{x,  ([>,  t)  =  ^  similarly  for  v  and  w),  Equation  (I.  16)  may 

be  rewritten  as 


0  0^ 

+  +  Uj^n(x,4>)Uij(x.<t>)4Wt)qij(t) 


where  mn  4  ij 

Since  the  generalized  coordinates  are  the  principal  or  normal  coordinates,  the  corre¬ 
sponding  vibrations  are  the  principal  modes  of  vibration,  and  the  products  of  the  velocities  in 
Equation  (1.  17)  vanish.  That  is,  from  Equations  (I.  12),  any  one  of  the  last  three  terms  in 
Equations  (I.  17)  would  be  of  the  form 


x=L 


q  q. .  r  F^(x)F-{x) 

^mn^ij  J  m'  '  i' 


x=0 


>  =  2^^ 
(j)  =  0 


dx 


where  Gj^(<t))  =  sin  n<j>  or  cos  n<t),  and  =  sin  j4>  or  cos  j(j>. 

Thus 


(b  =  Ztt 


G^(<t>)Gj((t>)dc|) 


(for  n  ^  j) 


and  Equation  (I.  17)  becomes 


T  =  Z  [qmn<t)]^[Umn  +  ^mn  + 


(I.  18) 


where 


-L 


"mn  =  /  I 


(1.19) 


0  0 


and  similarly  for  and 
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The  expression  for  the  potential  (strain)  energy  in  terms  of  the  derivative  of  the  dis¬ 
placement  components  and  the  generalized  coordinates  may  be  obtained  by  utilizing  Equations  (I.  2) 
and{I.  15)in  the  following 

Fa  ^+h/2  r  l-v^T 

V  =  —; - 27  I  II  Ux+  ^<|>+  2ve^e^  +  -7-  dxdcl,dz  (1.20 

2(1  -  )  -h/2  0  0  L  -* 

Since  terms  containing  products  of  the  generalized  coordinates  vanish,  Equation  (I.  20) 
may  be  rewritten  as  follows 

V  =  [•>// 

2(1  -  V  )  Loo 


.3  L  .Z-n  1 

+  j  J  G(u,  V,  w)dxd<j>  (1-21) 

0  0  J 

where  F(u,v,  w)  and  G(u,v,  w)  denote  a  collection  of  terms  involving  the  partial  derivatives  of 
u^^(x,  <()),  <|)),  and  <(>)  with  respect  to  x  and  <|>. 

From  Equation  (I.  21),  we  have 


Ztt 


F(u,  V,  w)dxd(j) 


av 


=  H(u,v,  w)q  (t) 


where  H(u,  v,  w)  are  merely  the  combined  constant  and  integral  terms. 
For  Lagrange's  equations,  we  have 


(1.22) 


_d  r  OT  1 


ax  av 

+ 


S'lmnd) 


Qmn(l) 


(1.23) 


From  Equation  (I.  20) 


ax 


=  0 


_d  r  ax  1 
dt  |^aq^,,(t)J 


p3.h<imn^l)I'*mn  ''^mn  ‘I  '*'mnl 


(I.  24) 


(I.  25) 


Introducing  Equations  (I.  22),  (I.  24),  and  (I.  25)  into  Equation  (I.  22) 


P^l^*mn  +  *mn  +  ^mJ  ^mn<l>  +  ^mn^D  " 


(I.  26) 
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Expressing  Equation  (I.  26)  in  a  more  convenient  form,  we  have 


where 


M 


{1.27) 


M  =  pah[u  +  V  +  ^TYinl 
mn  ^  *■  mn  mn  mn-" 


(generalized  mass) 


H(u,  V,  w) 
M _ 


Q _  (t)  =  Generalized  force 

mn 


L  2,r 

J  J  P(x,  (j),  t)  (x,  (j))  adxdtj) 
0  0 


For  the  numerical  solution  of  Equation  (I.  27),  the  method  of  finite  differences  was 
used,  and,  taking  the  second  central  differences,  we  have 


or 


qmn(ti) 


qmn<tr)  '  +  qmn{ti) 

(At)2 


‘lmn(tr) 


Q  (t.) 

mn'  i' 


(At) 2  +  [2  -  c.2^riAt)2]  q^n(ti)  -  qmn(tf) 


(1.28) 


where 


tr  =  ti  +  At 

t,  •=  t.  -  At 
i  1 

I.  3  Development  of  Load  Spectrums 


In  obtaining  the  expression  for  the  generalized  force  in  Equation  (I.  27),  use  was  made 
of  the  loading  information  for  cylindrical  shells  contained  in  References  I.  3  and  I.  4,  as  well  as 
other  publications.  Basically,  the  time  variation  of  the  loading  along  any  one  element  of  the 
cylinder  can  be  considered  of  the  form  shown  in  Figure  I.  3  or  I.  4,  the  former  utilizing  the  infor¬ 
mation  from  Reference  I.  3,  and  the  latter  from  Reference  I.  4. 

In  both  cases,  time  is  measured  from  the  instant  the  shock  front  strikes  the  leading 
edge  (element  at  (()  =  0®)  of  the  cylinder.  Again,  both  forms  include  the  arrival  time  (tg^)  for  each 
particular  cylinder  element,  the  duration  of  the  diffraction  loading  phase  (tj^)  and  the  positive 
phase  duration  (t.^.).  The  drag  phase  (for  t  >  tj^)  is  given  by 

P  =  Pi(t)  +  C^(<|>)q(?)  (1.29) 
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FIGURE  I.  3.  LOADING  SPECTRUM  AT  POINTS  ON 
FRONT  HALF  OF  CYLINDER 


FIGURE  I.  4.  SCHEMATIC  OF  PRESSURE  PROFILE  FOR  SHELL 
AT  ELEMENT  d>  =  (b. 


where,  for  the  time  increment  0  <  t  <  t+  measured  chronologically  from  t^,  the  relationship  for 
Pj  (t)  and  q(t)  is  given  in  Figure  I.  5,  and  C^i(4>)  is  the  drag  coefficient  for  the  particular  cylinder 
element  under  consideration.  For  the  loading  shown  in  Figure  I.  3,  the  drag  coefficient  Cjj(<|))  is 
obtained  from  Figure  6.  82b,  page  272,  Reference  I.  3;  for  the  loading  shown  in -Figure  I.  4,  C^jCcj)) 
is  found  from  the  curve  reproduced  in  Figure  I.  6. 

The  essential  difference  between  the  pressure  profiles  shown  in  Figures  I.  3  and  I.  4  lies 
in  the  diffraction  phase.  Since  the  pressure  profile  in  Figure  I,  4  represented  a  more  realistic 
loading  condition,  there  was  reason  to  tend  towards  this  approach.  Taking  into  acccount  the  experi¬ 
mental  model’s  dimensional  ratios  (L/D),  the  crossover  times  (D/U)  and  the  anticipated  (as  well  as 
experimentally  recorded)  durations,  there  appeared  to  be  adequate  aerodynamic  similitude  to 
justify  using  the  experimental  curves  in  Reference  I.  4  in  determining  the  appropriate  values  of 

tj  and  tj^  for  this  investigation's  loadings.  Moreover,  the  maximum  deformations  were  found 
to  occur  at  times  approximately  equal  to  or  less  than  the  diffraction  phase  duration  (tp).  This 
meant  that  the  diffraction  phase  constituted  the  major  (and,  in  some  cases,  the  entire)  deforma- 
tional  loading. 

The  curves  used  in  determining  the  significant  values  of  time  and  loading  required  to 
define  the  load  spectrum  in  Figure  I.  4  are  given  in  Figures  I.  7  through  1.9.  The  basic,  computa¬ 
tional  procedure  was  as  follows:  the  values  of  the  peak  incident  overpressure  (Pj)  and  positive 
phase  duration  (t^)  were  obtained  from  Tables  III.  1  and  III.  2.  The  variation  of  the  incident  over¬ 
pressure  and  the  dynamic  pressure  (see  Fig.  I.  5)  were  found  from  the  following  relations 


Pl(t^ 


q(t) 


=  Pi 


=  'Id 


2 


where 


‘Id 


%- 

2.  5Pj 

7p-r^ 


(psi) 


P^  =  Ambient  pressure 

Next,  the  crossover  time  (D/U)  was  computed,  where 


Cq  =  Ambient  sound  velocity  (ft  per  sec) 


(1.30) 


(1.31) 


(1.32) 


From  Figures  1.7  and  I.  8,  the  appropriate  values  for  t^,  t^  and  were  obtained  for  each  element 
in  the  cylinder’s  surface,  A  numerical  analysis  to  determine  the  optimum  number  of  elements  to 
be  considered  showed  that  for  increments  of  A<j)  =10°,  15°  or  22.  5°,  only  minor  differences 
existed  in  the  final  numerical  results.  Thus,  in  all  subsequent  computations,  the  elements  con¬ 
sidered  were  those  at  (j)  =  0°,  22.  5°,  45°,  .  .  .  ,  180°. 

The  values  of  Rj^((|>)  for  each  element  were  obtained  from  Figure  1.  9-  The  appropriate 
value  of  Pj  corresponding  to  the  time  t^  was  determined  from  Equation  (I.  30)  (where  t  =  t^)  or 
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from  Pi(t)  curves  such  as  are  shown  in  Figure  I.  5.  The  critical  pressures  were  then  found  for 
each  element 

P.  =  Ri(c|))Pj(tj) 

Typical  values  for  critical  pressures  and  times  are  given  in  Table  I.  1;  the  corresponding 
load  spectrums  are  shown  in  Figure  I.  10. 
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O  4  5  /3S  /<So 

OS/BSiTATlOKl  i  ,  C  0£a>/2.£B.S) 


FIGURE  I.  6.  VARIATION  OF  DRAG  COEFFICIENT  WITH  C SUNDER  ELEMENT 
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FIGURE  I.  7.  TIME  OF  ARRIVAL  OF  SHOCK  FRONT  ON  CYLINDER  SURFACE 


FIGURE  I.  8.  VARIATION  OF  CRITICAL  TIMES  (t^)  AND  DIFFRACTION  PHASE 
DURATION  (tjj)  WITH  CYLINDER  ELEMENT 
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ClVL/A/0£f2.  £L.aiA/f£/s/T  ( m 


FIGURE  I.  9.  VARIATION  OF  CRITICAL  PRESSURE  (P-)  TO  PEAK  INCIDENT  OVERPRESSURE  (Pj) 

RATIOS  WITH  CYLINDER  ELEMENT 


O  ./  .2  *3  •S  .7  .B  -9  /.e  /./  /.z  /.3  /.4  /,£ 


T-ZMe.  -  (^/99S£c.) 

FIGURE  I.  10.  PRESSURE  PROFILES  ALONG  VARIOUS  ELEMENTS  OF  CYLINDRICAL  SURFACE 
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TABLe4.  1.  CRITICAL  TIMES  AND  PRESSURES  (Pj  =  33  psi. 


=  1  msec)  D/U  =  0. 

53  msec 

msec 

^1 

msec 

Pi 

psi 

msec 

P2 

psi 

^3 

msec 

P3 

psi 

msec 

0 

0 

.  08 

85 

.  32 

18 

.  51 

21 

1.7 

22.  5 

.  02 

.  10 

70 

.  37 

14 

.  65 

9 

1.7 

45.  0 

.  07 

.  11 

55 

.41 

13 

.77 

3 

1.8 

67.  5 

.  16 

.  12 

39 

.  48 

11 

.97 

0 

2.  2 

90.  0 

.  27 

.  15 

24 

.  52 

10 

1.  1 

0 

3.0 

112.  5 

.  37 

.  16 

15 

.  58 

8 

1.3 

0 

3.  3 

135.0 

.  45 

.  19 

12 

.  62 

7 

1. 5 

0 

3.  3 

157.  5 

.  51 

.  20 

16 

.  68 

5 

1.6 

0 

3.  0 

180.0 

.53 

.  21 

25 

.74 

4 

1.8 

0 

2.7 
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APPENDIX  U 


EXPERIMENTAL  PROGRAM 


II.  1.  OBJECTIVES 


It  was  deemed  advisable  to  undertake  a  limited  experimental  program  which  would  provide 
data  to  test  the  adequacy  of  the  analytical  procedures.  Knowing  that  the  analytical  results  would  be 
in  the  form  of  predicted  displacements,  an  experimental  technique  was  evolved  for  obtaining  records 
of  the  response  of  cylindrical  shells  or  flat  plates  to  known  blast  loadings. 

II.  2.  MODEL  DESIGN 

II,  2.  1  Cylindrical  Shells 

All  the  cylindrical  shell  models  had  the  same  overall  nominal  dimensions:  diameter  -  12”, 
length  -  36”,  shell  wall  thickness  -  0.  036”.  The  three  basic  configurations  differed  with  respect  to 
the  end-closures  and  the  related  shell  boundary  condition.  For  the  Type  CA  specimen  (Figure  11.2), 
the  shell  was  welded  to  a  1/2”  thick  end  plate  so  as  to  closely  approximate  a  shell  with  fixed-end 
conditions.  The  end- closures  for  the  Type  CB  specimen  simulated  a  shell  with  a  simply- supported 
boundary  condition  (Figure  II.  3).  Type  CC  shell  models  use  end  plates  of  the  same  nominal  thick¬ 
ness  as  the  shell  (0.  036”).  These  plates  were  welded  to  the  shell  as  is  shown  in  Figure  II.  1.  The 
Type  CC  models  provided  an  elastic  boundary  and  were  used  to  demonstrate  the  effect  of  end 
restraints  of  flat  plate  end- closures. 

Within  each  CB  type  model,  a  scribe  plate  was  mounted  at  mid-length  on  the  three  support 
rods  (see  Figure  II.  3).  Four,  spring-loaded  styli,  supported  on  the  scribe  plate  and  bearing  against 
the  shell  wall,  were  used  to  record  the  deflection  at  mid- plane  points  on  the  0*,  90®,  180®,  and  270® 
elements  of  cylinder  (0®  denoting  the  blast  side  of  the  cylinder).  A  black  longitudinal  and  circumfer¬ 
ential  grid  on  one -inch  spacing  was  painted  on  the  yellow  surfaces  of  the  models  used  in  series  D,  A, 
G  and  H. 


II.  2.  2.  Circular  Flat  Plates 


The  test  fixture  used  to  support  the  circular  flat  plates  is  shown  in  Figure  11.4.  As  shown 
in  the  drawing,  the  plate  is  simply  supported  at  the  boundary;  for  the  fixed  boundary  condition,  an 
oversized  plate  was  used  and  clamped  between  the  inner  pipe  and  the  outer  ring.  In  order  to  evaluate 
the  effects  of  membrane  and  beam  action,  0.  0312”  and  0.  0598”  thick  plates  were  tested.  White 
radial  lines  at  45®  were  painted  on  blue  plate  surfaces. 

II.  3.  INSTRUMENTATION 


Instrumentation  consisted  of  synchronized,  high-speed  photography  to  record  the  cylinders' 
dynamic  response,  mechanical  styli  recordings  of  maximum  deformations,  and  a  quartz  transducer 
system  to  record  overpressures  and  positive  phase  durations.  Each  model  was  covered  by  one  or 
two  high-speed  cameras  with  film  speeds  from  3,  000  to  5,  000  frames  per  second.  A  sequence  cir¬ 
cuit  was  constructed  to  synchronize  the  detonation  and  camera  action  so  that  the  initial  impact  of  the 
shock  wave  on  the  model  would  occur  when  the  cameras  had  obtained  their  maximum  speed. 

Two  quartz  transducers  were  flush- mounted  in  a  18"  X  24”  steel  surface  plate.  An  adap¬ 
tor,  consisting  of  a  double  brass  fitting  separated  by  a  viscous  potting  material,  was  constructed  to 
minimize  ground  acceleration  effects  in  the  pickup.  The  transducers  were  placed  at  ground  level  at 
radial  distances  from  ground  zero  equal  to  those  from  ground  zero  to  the  model.  The  transducers' 
outputs  were  fed  through  arhplifier- calibrators  to  an  oscilloscope  and  recorded  by  oscillograph 
cameras . 

n.4.  TEST  PROGRAM 


A  total  of  nineteen  charges  in  nine  series  was  detonated.  The  charges  used  and  the  cylin 
drical  shell  models  and  flat  plates  exposed  in  each  shot  are  given  in  Table  II.  1.  Each  cylindrical 
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TABLE  II.  1.  SUMMARY  OF  EXPERIMENTAL  PROGRAM 
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P2  --  5  0.25  Not  useable 


shell  model  was  suspended  with  nylon  cords;  the  position  and  orientation  of  each  model  relative  to 
ground  zero  are  shown  in  Figures  II.  5  through  II.  8;  the  controlling  distances  are  given  in  Table  II.  1. 

Shot  series  A  through  C  were  used  to  obtain  qualitative  information  on  the  deformational 
response  of  the  cylindrical  shells  as  influenced  by  such  factors  as  the  distance  and  orientation  of 
models  with  respect  to  ground  zero.  In  particular,  the  tests  were  used  to  determine  the  modes  of 
failure  (in  the  form  of  permanent  deformation  or  fracture)  as  related  to  estimated  overpressure. 

The  tests  using  5  lb  of  flaked  TNT  (Series  E  and  G)  were  conducted  on  Southwest  Research 
Institute  grounds  for  the  purpose  of  checking  out  the  equipment  and  instrumentation.  In  Series  G,  the 
distance  between  the  model  and  the  charge  was  10  feet,  and,  at  this  distance,  it  was  possible  to 
photograph  small  (but  measurable)  elastic  deformations  and  mode  shapes.  The  ground  zero-to- 
model  distances  for  the  30-lb  shots  at  Camp  Bullis  varied  from  Z6  feet  to  15  feet;  each  series  of 
shots  terminated  with  failure  of  the  model  in  the  form  of  large,  permanent  deformations.  (See 
Figures  11.  9  through  II.  1  7 . ) 
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FIGURE  II.  5.  MODEL  POSITIONS  FOR  SHOTS  Al,  A2.  AND  A3 


FIGURE  II.  6.  MODEL  POSITIONS  FOR  SHOTS  Bl,  B2,  AND  B3 


FIGURE  II.  9.  END  PLATE  FAILURE  IN 
MODEL  CCl  AFTER  SHOT  A3 
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FIGURE  II.  10.  END  PLATE  AND  SHELL  FAILURE 
IN  MODEL  CCl  AFTER  SHOT  B3 
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FIGURE  II.  1 1.  SLIGHT  REDUCTION  IN  SHELL  DIAMETER 
IN  MODEL  CC2  AFTER  SHOT  BI 


FIGURE  II.  12.  SHELL  AND  END  PLATES  PERMANENT 
DEFORMATIONS  IN  MODEL  CC2  AFTER  SHOT  B2 


FIGURE  II.  13.  SHELL  AND  END  PLATE  PERMANENT 
DEFORMATIONS  IN  MODEL  CC3  AFTER  SHOT  Cl 


FIGURE  n.  14.  END  PLATE  DEFORMATIONS  IN  MODEL  CC4 
AFTER  SHOT  C7 


Vi 

0 

FIGURE  11.  15.  END  PLATE  DEFORMATION  IN  MODEL  CC5 
AFTER  SHOT  8 


FIGURE  II.  16.  SHELL  DEFORMATION  IN  MODEL  CA2 
AFTER  SHOT  7 
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APPENDIX  m 


EXPERIMENTAL  RESULTS 
in.  1 .  PRESSURE- TIME  MEASUREMENTS 


In  developing  the  experimental  portion  of  the  program,  consideration  was  given  to  creating 
blast  phenomena  with  conventional,  chemical  high  explosives  which,  when  applied  in  an  analysis  pro¬ 
cedure  to  scaled  models  of  missile  bodies,  would  provide  a  discernible  equivalence  for  the  nuclear 
weapon  effects  on  actual  missile  structures.  To  this  end,  full  use  was  made  of  the  information  avail¬ 
able  in  Reference  III.  1,  as  well  as  other  supporting  publications  containing  information  on  HE  explo¬ 
sions  . 


The  peak  shock  overpressure,  as  related  to  distance,  plays  a  leading  role  in  explosion- 
damage  correlations.  Selection  of  the  controlling  parameters  for  the  HE  blast  phenomena,  there¬ 
fore,  began  with  the  relation  shown  in  Figure  III.  1  between  the  peak  overpressure  in  free  air  and 
scaled  (radial)  range  as  proposed  by  the  Kirkwood- Brinkley  theory(^^' and  substantiated  by  the 
experimental  evidence  of  Fisher(^^^‘ and  Weibull^^^^' In  this  relation,  the  scaled  range  is 
defined  as 


(W)l/3 


where 


R  is  the  radius  (feet) 

W  is  the  weight  (yield)  (pounds) 

Since  the  HE  explosions  would  take  place  on  or  near  the  ground,  other  blast  phenomena 
associated  with  surface  or  near  surface  bursts  were  taken  into  account.  Foremost  among  these 
considerations  was  the  ground  reflection  of  the  blast  wave  and  the  subsequent  coalescence  of  the 
incident  and  reflected  waves  in  the  development  of  the  Mach  stem. 


In  accordance  with  one  definition,  let  be  the  size  of  burst  exploded  over  a  reflector  at 
a  height  of  burst  equal  to  hp.  At  a  slant  range  of  Rp  >  2  hp  (thereby  placing  the  point  in  question  in 
the  faV  Mach  region),  the  reflected  peak  overpressure  is  Pp.  For  a  free -air  burst  of  Wj,  the  inci¬ 
dent  (or  free -air)  peak  overpressure  at  a  range  of  Rj  equal  to  Rp  would  be  Pj.  Let  the  yield  of  Wj 
be  such  that,  at  a  range  of  Rj  =  Rp,  the  incident  and  reflected  peak  overpressures  are  equal.  Since 
the  scaled  slant  ranges  are 


Xp  = 


(Wr) 


1/3 


% 

(Wj)l/3 


by  definition 

Wi 

Rp  (reflection factor)  =  — — 

Wr 


3 


Accordingly,  for  a  given  Wp  at  height  of  burst  hp,  one  would  expect  at  a  range  of  Rp  a  reflected 
peak  overpressure  in  the  Mach  reflection  region  equal  to  that  obtained  from  a  free -air  burst  with 
a  yield  of  RpWj. 

An  indication  of  the  reflection  factor  to  be  expected  in  the  Mach  reflection  was  obtained 
from  the  information  presented  in  Figure  III.  2  relating  the  reflection  factor  for  TNT  to  the  scaled 
height  of  burst  for  a  particular  (and  not  necessarily  ideal)  reflecting  surface.  This  and  the  additional 
information  for  pentolite  and  the  related  data  for  nuclear  bursts  indicated  that  scaled  heights  of 
bursts  could  be  selected  such  that  the  reflection  factor  in  the  far  Mach  region  would  be  between  0.,5 
and  2.0,  approximately.  In  Figure  III.  1  are  shown  the  R  W  curves  for  these  two  extreme s  relating 
the  reflected  peak  overpressure  with  the  scaled  slant  range.  The  curves  in  Figure  III.  3  provide  the 
same  information  for  a  nonscale d  range. 
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FIGURE  III.  1.  PEAK  OVERPRESSURE  VS  SCALED  DISTANCE  FOR 
SPHERICAL  TNT  IN  FREE -AIR  AT  SEA  LEVEL 
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FIGURE  III.  2.  REFLECTION  FACTORS  FOR  TNT  EXPLOSIONS 
AT  LOW  HOB  OVER  VARIOUS  SURFACES 
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FIGURE  III.  3.  PEAK  OVERPRESSURE  VERSUS  RANGE  FOR  TNT 
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Since  the  above  method  of  predicting  the  reflected  peak  overpressure  was  limited  to  the  far 
Mach  region  and  in  order  to  avoid  having  the  missile  model  experience  two  pressure  peaks  in  the 
regular  reflection  region,  the  charge  size,  height  of  burst  and  the  ground  range  were  selected  such 
that  the  triple  point  would  be  above  the  uppermost  portion  of  any  model,  regardless  of  its  orientation. 

The  information  in  the  literature  provided  the  necessary  relation  between  the  Mach  stem 
height  and  the  range  for  various  heights  of  burst.  From  prior  qualitative  tests,  it  was  known  that  the 
missile  models  would  sustain  damage  from  a  30-lb  charge  of  TNT  at  a  ground  range  of  1  5  to  Z6  feet, 
depending  on  the  orientation  and  construction  of  the  model  itself.  The  model  positioned  with  its  lon¬ 
gitudinal  axis  perpendicular  to  ground  represented  the  most  severe,  triple-point  height  requirement. 
For  this  arrangement  and  a  ground  range  of  1  2  feet,  the  scaled  height  of  burst  would  be  limited  to  a 
scaled  HOB  of  0.  8  or  less  or,  for  a  charge  weight  of  30  lb  of  TNT,  a  burst  elevation  above  ground  of 
2-1/2  feet  or  less.  Referring  to  Figure  III.  2,  this  height  of  burst  would  provide  a  reflection  factor 
of  approximately  1.8  and  therefore  ground  reflection,  peak  overpressures  at  the  model  of  from 
approximately  15  to  100  psi  depending  on  the  ground  range. 

It  was  recognized  that  such  factors  as  the  actual  ground  leflection  characteristics  at  the  test 
site,  the  type  of  charge  and  the  charge  shape  would  all  somewhat  influence  the  blast  phenomena. 
Accordingly,  pressure  measuring  devices  were  used  to  determine,  among  other  things,  the  extent 
of  the  deviation  (if  any)  of  the  predicted  peak  overpressure.  (See  Appendix  II.  ) 

The  results  of  these  pressure  readings  (shown  as  points  superimposed  on  the  curves  shown 
in  Figures  III.  1  andlll.  2)  are  especially  encouraging  particularly  in  view  of  the  uncertainties 
normally  associated  with  blast  phenomena  measurements.  It  is  seen  that  with  the  exception  of  the 
30-lb,  H4  shot  and  the  5-lb,  El  and;E2  shots,  a  ground  reflection  factor  of  between  1. 6  and  1.  7  for 
the  30-lb  shots  and  of  0.  7  for  the  5-lb  shots  is  quite  valid.  The  values  for  Rp  were  used  for  predict¬ 
ing  the  peak  reflected  overpressures  in  those  tests  where  pressure  measurements  were  not  obtained. 
Table  III.  1  summarizes  the  pertinent  peak  overpressure  and  reflection  factor  data. 

In  addition  to  the  overpressures,  the  other  blast  parameters  which  significantly  influence  the 
extent  to  which  the  model  is  damaged  are  the  positive  phase  duration  and  the  positive  impulse.  The 
time  of  arrival  of  the  blast  wave  is  of  importance  in  the  establishment  and  interpretation  of  the 
model's  high-speed  photography  data. 

^Available  free -air  curves  for  arrival  time,  positive  phase  duration  and  positive  impulse  for 
TNT  provided  suitable  approximations  for  adjusting  the  pressure  recording  and  high-speed  photo¬ 
graphy  equipment.  Table  III.  2  offers  a  comparison  between  measured  and  predicted  data.  A  sum¬ 
mary  of  the  pressure  data  is  given  in  Table  III.  3. 

III.  2.  HIGH-SPEED  PHOTOGRAPHY  DATA 


The  variation  of  the  radial  displacement  at  x  =  L/2,  (|)  =  0®  with  time  (where  t  =  0  denotes  the 
arrival  of  the  blast  wave  at  the  cylinder's  leading  edge)  is  shown  in  Figures  III.  4  through  III.  13. 
Each  of  the  curves  was  fitted  to  points  obtained  from  either  prints  or  projections  of  each  individual 
frame . 


The  curves  in  Figures  III.  8  and  III.  10  need  some  additional  interpretation.  The  model  (CB2) 
used  in  shot  HI  had  been  previously  subjected  to  four,  5-lb  explosions  in  shot  series  G.  These 
smaller  (but  repeated  exposures)  resulted  in  a  slight  permanent  set  at  the  cylinder's  top,  leading 
edge  adjacent  to  the  grooved  (for  simple  support)  end  plates.  In  the  subsequent  shot  HI,  this  edge 
of  the  shell  moved  out  of  the  groove  at  approximately  1. 7  msec  (see  Figure  III.  8)  after  the  arrival 
of  the  blast  wave.  In  so  doing,  the  subsequent  response  of  the  shell  was  no  longer  that  of  a  cylinder 
with  simple  supports  at  each  boundary.  However,  prior  to  that  time,  the  shell's  response  followed 
the  theoretically  predicted  displacement  pattern. 

The  same  mode  of  failure  was  experienced  by  Model  CB3  (Figure  III.  10)  in  shot  H4.  Here, 
the  shell's  leading  edge  at  the  top  end- plate  began  to  move  out  of  the  groove  at  approximately  1.  5 
msec.  Although  the  shell  did  subsequently  return  to  zero  displacement  at  x  =  L/2,  cj)  =  0®,  the  condi¬ 
tion  of  simple  support  was  no  longer  in  effect.  For  this  reason,  the  simple- support  response  was 
limited  to  that  portion  of  the  displacement  curve  prior  to  1 . 5  msec. 
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Shot 

No. 

El 

KZ 

E3 

FI 

FZ 

G1 

G2 

G3 

G4 

H2 

H3 

H4 

J1 

JZ 

J3 


N.  A. 

(1) 

(^) 

(3) 

(A) 


(B) 

(C) 

(D) 

(E) 


TABLE  III.  1.  PEAK  OVERPRESSURE  IN  MACH  REFLECTION 
REGION  (EXPERIMENTAL  DATA) 


Ground^ 


Charge^  ^ ) 
(lb) 

hob(^) 

(ft) 

Range 

(ft) 

(A) 

RF  =  0.  5 

rf  = 

2.  0 

(B) 

Rf 

(C) 

(D) 

(E) 

5 

0. 

33 

10 

12.  8 

33. 

0 

1. 

32 

25. 

1 

>5. 

0 

2.  0 

5 

0. 

33 

10 

12.  8 

33. 

0 

1. 

32 

25. 

1 

29. 

3 

1.  6 

5 

0. 

33 

10 

12.  8 

33. 

0 

1. 

32 

25. 

1 

18. 

8 

0.9 

30 

1. 

5 

20 

10.  5 

27. 

2 

1. 

89 

27. 

U 

23. 

0 

1.  6 

30 

1, 

5 

16 

16.  3 

47. 

0 

1. 

89 

46. 

4 

41. 

0 

1.  7 

5 

0. 

33 

10 

12.  8 

33. 

0 

1. 

32 

25. 

1 

19. 

7 

0.9 

5 

0. 

33 

10 

12.  8 

33. 

0 

1. 

32 

25. 

1 

16. 

1 

0.  7 

5 

0. 

33 

10 

12.  8 

33. 

0 

1. 

32 

25. 

1 

16. 

5 

0.  7 

5 

0. 

33 

10 

12.  8 

33. 

0 

1. 

32 

25. 

1 

N.  i 

N.  A. 

30 

2. 

5 

20 

10.  5 

27. 

2 

1. 

89 

27. 

0 

24. 

6 

0.  7 

30 

2. 

5 

20 

10.  5 

27. 

2 

1. 

89 

27. 

0 

23. 

0 

1.  6 

30 

2. 

5 

15 

18.  5 

55. 

9 

1. 

89 

54. 

0 

53. 

0 

1.  6 

5 

0. 

33 

10 

12.  8 

33. 

0 

1. 

32 

25. 

1 

13. 

1 

0.  6 

5 

0. 

33 

10 

12.  8 

33. 

0 

1, 

32 

25. 

1 

16. 

2 

0.  7 

5 

0. 

33 

5 

39.  2 

140. 

0 

1. 

32 

102. 

0 

64. 

0 

0.  7 

Not  available. 

Flaked  TNT. 

Height  of  burst. 

Distance  from  ground  zero  to  model  and/or  pressure  pickup. 

Estimated  range  of  peak  overpressure  in  Mach  reflection  region  (psi); 
Kirkwood-Brinkley  theory  -  see  Figures  III.  1  and  III.  3. 

Reflection  factor  for  hard  clay  surface  as  function  of  HOB  (see  Figure  III.  2). 
Peak  overpressure  based  on  Rp  from  Column  (B). 

Experimentally  measured  peak  overpressure  (psi). 

Adjusted  reflection  factor  based  on  measured  peak  overpressure. 
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TABLE  m.  2.  TIME  OF  ARRIVAL.  POSITIVE  PHASE  DURATION 
IN  MACH  REFLECTION  REGION  (EXPERIMENTAL  DATA) 


Shot 

Charge^^) 

Ground^^) 

(A) 

(B) 

(C) 

(D) 

No. 

(lb) 

(ft) 

Range  (ft) 

(msec) 

(msec) 

(msec) 

(msec) 

D1 

30 

2.5 

26 

11.5 

N.A. 

_  . 

_  -  - 

D2 

30 

2.5 

22 

9.3 

7.9 

— 

--- 

D3 

30 

2.5 

16 

5.  0 

5.3 

— 

— 

El 

5 

0.  25 

10 

3.4 

N.A. 

2.  6 

2.  7 

E2 

5 

0.  25 

10 

3.4 

N.A. 

2.6 

2.  8 

E3 

5 

0.  25 

10 

3.4 

N.A. 

2.  6 

2.2 

FI 

30 

1.5 

20 

7.4 

N.A. 

5.  1 

4.2 

F2 

30 

1.5 

16 

4.9 

N.A. 

4.  0 

4.3 

G1 

5 

0.  25 

10 

3.4 

4.  5 

2.  6 

2.4 

G2 

5 

0.25 

10 

3.4 

3.9 

2.  6 

2.  3 

G3 

5 

0.25 

10 

3.4 

3.9 

2.  6 

2.  5 

G4 

5 

0.25 

10 

3.4 

3.8 

2.  6 

2.  7 

HI 

30 

2.  5 

20 

7.4 

7.  1 

5.  1 

N.A. 

H2 

30 

2.  5 

20 

7.4 

N.A. 

5.  1 

6.  0 

H3 

30 

2.  5 

20 

7.4 

N.A. 

5.  1 

«4.  5 

H4 

30 

2.  5 

15 

3.  7 

N.A. 

3.9 

N.A. 

J1 

5 

0.  25 

10 

3.4 

N.A. 

2.  6 

3.  1 

J2 

5 

0.  25 

10 

3.4 

N.A. 

2.  6 

3.4 

J3 

5 

0.  25 

5 

1.2 

N.A. 

1.6 

1.  1 

N.A.  Not  available. 

(1)  Flaked  TNT. 

(2)  Height  of  burst. 

(3)  Distance  from  ground  zero  to  model  and/or  pressure  pickup. 

(A)  Estimated  time  of  arrival  -  tg. 

(B)  Recorded  time  of  arrival  -  t^  (from  high  speed  photography). 

(C)  Estimated  positive  phase  duration  -  t+ 

(D)  Recorded  positive  phase  duration  -  t-\. 


TABLE  m.  3.  SUMMARY  OF  PEAK  INCIDENT  OVERPRESSURES 
AND  POSITIVE  PHASE  DURATIONS 


(1) 

E:  end- on  loading;  S:  side- on 

(4) 

E  -  elastic;  F  -  fixed 

loading;  O:  oblique  loading 

S  -  simple  support 

(2) 

Peak  incident  overpressure 

(5) 

Experimental  -  See 
Table  m.  1  or  III.  2 

(3) 

Positive  phase  duration 

Shot 

Model  (1) 

Pi  (2) 

t+  (3) 

Boundary- 

Model 

No. 

Orientation 

(psi) 

(msec) 

Restraint  (4) 

CCl 

A1 

E 

4 

7 

E 

A2 

10 

6 

A3 

24 

4 

B1 

c 

) 

18 

5 

B2 

23 

4.  5 

1 

B3 

32 

4 

CC2 

A1 

2 

9 

A2 

6 

7 

A3 

17 

5 

B1 

29 

4 

B2 

38 

3 

CC3 

Cl 

31 

5 

CC4 

Cl 

I 

r 

22 

6 

CC4 

C2 

E 

27 

5 

CCS 

Cl 

O 

13 

7 

CCS 

C2 

O 

17 

6 

CAl 

Cl 

E 

22 

6 

1 

r 

CAl 

C2 

E 

27 

4 

CA2 

Cl 

£ 

27 

4 

CA3 

D1 

15 

6.5 

D2 

19 

5.  5 

D3 

42 

3.8 

CBi 

FI 

27 

4 

S 

CBl 

F2 

46 

4 

CB2 

G1 

20  (5) 

2.4  (5) 

G2 

16  (5) 

2.3  (5) 

G3 

17  (5) 

2.3  (5) 

G4 

18 

2.  5 

HI 

25 

5.  1 

CB3 

H3 

23  (S) 

4.5 

CB3 

H4 

53  (5) 

3.9 

PI 

J1 

- 

25 

3 

P2 

J2 

- 

25 

3 

P2 

J3 

- 

102 
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C> /■:>/:>/.  A  c£^/yfS'/^r  -  <rsrj  , 

SAO/AL  O/sfi^CACS/^e-f^T  -  y/  Q/Ae:./><S 


Tt/^AE. 


FIGURE  III.  4.  RADIAL  DISPLACEMENTS  AT  x  =  L/Z,  <t>  =  0  “  OBTAINED  FROM 
HIGH  SPEED  PHOTOGRAPH  OF  CYLINDRICAL  MODELS 


■Z  3 

-r/AA£.  (^AOS<=.c. ) 


FIGURE  (II.  5.  RADIAL  DISPLACEMENTS  AT  x  =  L/2,  <()  =  0 ’’  OBTAINED  FROM 
HIGH  SPEED  PHOTOGRAPH  OF  CYLINDRICAL  MODELS 
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/^AO/AL  /:>f  :sp^Ac:SM£^r 


FIGURE  HI  6  RADIAL  DISPLACEMENTS  AT  x  =  L/2,  <()  ^  0  °  OBTAINED  FROM 
HIGH  SPEED  PHOTOGRAPH  Of'  CYLINDRICAL  MODELS 


FIGURE  ITT,  7.  RADIAL  DISPLACEMENTS  AT  x  :=  L/2,  <()  =  0°  OBTAINED  FROM  HIGH 
SPEED  PHOTOGRAPH  OF  CYLINDRICAL  MODELS 


98 


/2/iO/AL  O/ 


2  3  B  ^  7 

T^/A^S  -  S<r<= ) 

FIGURE  m.  8.  RADIAL  DISPLACEMENTS  AT  x  =  L/2,  <|)  =  0”  OBTAINED  FROM 
HIGH  SPEED  PHOTOGRAPH  OF  CYLINDRICAL  MODELS 


/  2  3  ^  B  <»  7 


T/A/fB  - 

FIGURE  III.  9.  RADAL  DISPLACEMENTS  AT  x  =  L/2,  <})  =  0  “  OBTAINED  FROM 
HIGH  SPEED  PHOTOGRAPH  OF  CYLINDRICAL  MODELS 
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/:?/ S/=’L^C£‘/\>f£A/T 


CZ^/S  2 


MGO€L 
SJ^OT 

=  S  3  p  ^  /90  S^  CL^ 

(2)  To^  a^  czy^/A/z>£^ 
oc/r  £>/^  B,/yo~ /^i.Ar£ 


'T‘/Ayf£  - 

FIGURE  III.  10.  RADIAL  DISPLACEMENTS  AT  x  =  L/2,  <|)  =  0”  OBTAINED  F^OM 
HIGH  SPEED  PHOTOGRAPH  OF  CYLINDRICAL  MODELS 


MOa^L  CZB-l 
SHOT  O/ 


■,4- 


<S-3  r?7SOO 


T/S^£,  -  SO<^ 

FIGURE  m.  11.  RADIAL  DISPLACEMENTS  AT  x  =  L/2.  <|)  =  0*  OBTAINED  FROM 
HIGH  SPEED  PHOTOGRAPH  OF  CYLINDRICAL  MODELS 


[w] 


i 


k 

I 

VI 


5 

0 


3  Cc  7 

FIGURE  III.  IZ.  RADIAL  DISPLACEMENTS  ATx=L/Z,  <j)  =  0°  OBTAINED  FROM 
HIGH  SPEED  PHOTOGRAPH  OF  CYLINDRICAL  MODELS 


<*1 

VI' 

V) 


k 

? 

11 

0 

X 

>1 


vl 


5^0  T  D  3 
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/ 

\ 

■ 

■ 

/ 

\ 

\ 

ll 

n 

■ 

■ 

■ 

■ 

■ 

■ 

■ 

■ 

O  /  2  3  4-  s  y 


T/K^B  -  (^/n£:£.c:.^ 

FIGURE  ni.  13.  RADIAL  DISPLACEMENTS  AT  x  =  L/2,  <})  =  0°  OBTAINED  FROM 
HIGH  SPEED  PHOTOGRAPH  OF  CYLINDRICAL  MODELS 
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